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MATRIX-VALUED ORTHOGONAL POLYNOMIALS RELATED TO THE 
QUANTUM ANALOGUE OF (SU(2) x SU(2),diag) 

NOUD ALDENHOVEN, ERIK KOELINK, PABLO ROMAn 


Abstract. Matrix-valued spherical functions related to the quantum symmetric pair for the quantum 
analogue of (SU(2) x 311(2), diag) are introduced and studied in detail. The quantum symmetric pair is 
given in terms of a quantised universal enveloping algebra with a coideal subalgebra. The matrix-valued 
spherical functions give rise to matrix-valued orthogonal polynomials, which are matrix-valued analogues 
of a subfamily of Askey-Wilson polynomials. For these matrix-valued orthogonal polynomials a number 
of properties are derived using this quantum group interpretation: the orthogonality relations from the 
Schur orthogonality relations, the three-term recurrence relation and the structure of the weight matrix 
in terms of Chebyshev polynomials from tensor product decompositions, the matrix-valued Askey- 
Wilson type ^-difference operators from the action of the Casimir elements. A more analytic study of 
the weight gives an explicit LDU-decomposition in terms of continuous g-ultraspherical polynomials. 
The LDU-decomposition gives the possibility to find explicit expressions of the matrix entries of the 
matrix-valued orthogonal polynomials in terms of continuous g-ultraspherical polynomials and g-Racah 
polynomials. 
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1. Introduction 

Shortly after the introduction of quantum groups it was realised that many special functions of basic 
hypergeometric type [16] have a natural relation to quantum groups, see e.g. [9] Ch. 6], [21], [31] for 
references. In particular, many orthogonal polynomials in the g-analogue of the Askey scheme, see 
e.g. [26] ■ have found an interpretation on compact quantum groups analogous to the interpretation 
of orthogonal polynomials of hypergeometric type from the Askey scheme on compact Lie groups and 
related structures, see e.g. m. sa¬ 
in case of the harmonic analysis on classical Gelfand pairs one studies spherical functions and related 
Fourier transforms, see m- For our purposes a Gelfand pair consists of a Lie group G and a compact 
subgroup K so that the trivial representation of K in the decomposition of any irreducible representation 
of G restricted to K occurs with multiplicity at most one. The spherical functions are functions on G 
which are left- and right-A-invariant. The zonal spherical functions are realised as matrix elements of 
irreducible G-representations with respect to a fixed A-vector. For special cases the zonal spherical 
functions can be identified with explicit special functions of hypergeometric type, see m Ch. 9], m 
§IV]. The zonal spherical functions are eigenfunctions to an algebra of differential operators, which 
includes the differential operator arising from the Gasimir operator in case G is a reductive group. 
For special cases with G compact, we obtain orthogonality relations and differential operators for the 


Date-. July 14, 2015. 


1 








spherical functions, which can be identified with orthogonal polynomials from the Askey-scheme. For 
the special case G = SU(2) x SU(2) with K = SU(2) embedded as the diagonal subgroup, the zonal 
spherical functions are the characters of SU(2), which are identified with the Chebyshev polynomials 
Un of the second kind by the Weyl character formula. The quantum analogue of the Gelfand pair 
situation has been generalised to the setting of quantum groups, mainly in the compact context, see e.g. 
Andruskiewitch and Natale [3] for the case of finite-dimensional Hopf-algebra with a Hopf sub-algebra, 
Floris [a, Koornwinder [32) . Vainermann |44] for more general compact quantum groups, and, for a 
non-compact example, Caspers [B]. 

The notions of matrix-valued and vector-valued spherical functions have already emerged at the be¬ 
ginning of the development of the theory of spherical functions, see e.g. |15j and references given there. 
However, the focus on the relation with matrix-valued or vector-valued special functions only came later, 
see e.g. references given in [TB], [ID]. Griinbaum, Pacharoni, Tirao m give a group theoretic approach 
to matrix-valued orthogonal polynomials emphasising the role of the matrix-valued differential operators, 
which are manipulated in great detail. The paper [T7] deals with the case (G, K) = (SU(3), U(2)). Moti¬ 
vated by m and the approach of Koornwinder |29] , the group theoretic interpretation of matrix-valued 
orthogonal polynomials on (G, K) = (SU(2) x SU(2), SU(2)) is studied from a different point of view, in 
particular with less manipulation of the matrix-valued differential operators, in |23j . [24) . see also m. 
[40) . The point of view is to construct the matrix-valued orthogonal polynomials using matrix-valued 
spherical functions, and next using this group theoretic interpretation to obtain properties of the matrix¬ 
valued orthogonal polynomials. This approach for the case {G,K) = (SU(2) x SU(2),SU(2)) leads to 
matrix-valued orthogonal polynomials for arbitrary size, which can be considered as analogues of the 
Ghebyshev polynomials of the second kind. A combination of the group theoretic approach and analytic 
considerations then allows us to understand these matrix-valued orthogonal polynomials completely, i.e. 
we have explicit orthogonality relations, three-term recurrence relations, matrix-valued differential oper¬ 
ators having the matrix-valued orthogonal polynomials as eigenfunctions, expression in terms of Tirao’s 
[42) matrix-valued hypergeometric functions, expression in terms of well-known scalar-valued orthog¬ 
onal polynomials from the Askey-scheme, etc. This has been analytically extended to arbitrary size 
matrix-valued orthogonal Gegenbauer polynomials [33], see also [31] for related 2 x 2-cases. 

The interpretation on quantum groups and related structures leads to many new results for special 
functions of basic hypergeometric type. In this paper we use quantum groups in order to obtain matrix¬ 
valued orthogonal polynomials as analogues of a subclass of the Askey-Wilson polynomials. In particular, 
we consider the Ghebyshev polynomials of the second kind, recalled in (15.61) . as a special case of the Askey- 
Wilson polynomials |U (2.18)]. Moreover, we know that the Ghebyshev polynomials occur as characters 
on the quantum SU(2) group, see H?] §A.l]. The approach in this paper is to establish the quantum 
analogue of the group-theoretic approach as presented in m, m, see also m, m, for the example of 
the Gelfand pair G = SU(2) x SU(2) with K = SU(2). For this approach we need Letzter’s approach 
[34] . [35] . [36] to quantum symmetric spaces using coideal subalgebras. We stick to the conventions as 
in Kolb [27] and we refer to [281 §1] for a broader perspective on quantum symmetric pairs. So we 
work with the quantised universal enveloping algebra Uq{Q) = Uq{su{2) © su(2)), introduced in Section 
[31 equipped with a right coideal subalgebra B, see Section [3] Once we have this setting established, the 
branching rules of the representations of Uq{Q) restricted to B follow by identifying B with the image 
of Uq{5u(2)) (up to an isomorphism) under the comultiplication using the standard Clebsch-Gordan 
decomposition. In particular, it gives explicit intertwiners. Next we introduce matrix-valued spherical 
functions in Section |3J Using the matrix-valued spherical functions we introduce the matrix-valued 
orthogonal polynomials. Then we use a mix of quantum group theoretic and analytic approaches to 
study these matrix-valued orthogonal polynomials. So we find the orthogonality for the matrix-valued 
orthogonal polynomials from the Schur orthogonality relations, the three-term recurrence relation follows 
from tensor-product decompositions of 77q(g)-representations, the matrix-valued g-difference operators 
for which these matrix-valued orthogonal polynomials are eigenvectors follow from the study of the 
Gasimir elements in lAq{Q). More analytic properties follow from the LDU-decomposition of the matrix¬ 
valued weight function, and this allows to decouple the matrix-valued g-difference operators involved. 
The decoupling gives the possibility to link the entries of the matrix-valued orthogonal polynomials 
with (scalar-valued) orthogonal polynomials from the g-analogue of the Askey-scheme, in particular the 
continuous q-ultraspherical polynomials and the q-Racah polynomials. The approach of m does not 
seem to work in the quantum case, because the possibilities to transform g-difference equations are very 
limited compared to transforming differential equations. We note that in [3] §5] matrix-valued spherical 
functions are considered for finite-dimensional Hopf algebras with respect to a Hopf-subalgebra. 
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The approach to matrix-valued orthogonal polynomials from this quantum group setting also leads 
to identities in the quantised function algebra. This paper does not include the resulting identities 
after using infinite-dimensional representations of the quantised function algebra. Furthermore, we have 
not supplied a proof of Lemma 15.41 using infinite-dimensional representations and the direct integral 
decomposition of the Haar functional, but this should be possible as well. 

In general, the notion of a quantum symmetric pair seems to be best-suited for the development of 
harmonic analysis in general and in particular of matrix-valued spherical functions on quantum groups, 
see e.g. EH], ED, ESI, ESI, EZI, EH] and references given there. When considering other quantum 
symmetric pairs in relation to matrix-valued spherical functions, the branching rule of a representation 
of the quantised universal enveloping algebra to a coideal subalgebra seems to be difficult. In this paper 
it is reduced to the Clebsch-Gordan decomposition, and there is a nice result by Oblomkov and Stokman 
[381 Prop. 1.15] on a special case of the branching rule for quantum symmetric pair of type AIII, but in 
general the lack of the branching rule for the quantum symmetric pairs is an obstacle for the study of 
quantum analogues of matrix-valued spherical functions of e.g. m, EZ], m, EH]- 

The matrix-valued orthogonal polynomials resulting from the study in this paper are matrix-valued 
analogues of the Chebyshev polynomials of the second kind viewed as an example of the Askey-Wilson 
polynomials. We expect that it is possible to obtain matrix-valued analogues of the continuous q- 
ultraspherical polynomials viewed as subfamily of the Askey-Wilson polynomials using the approach 
of ES] using the Askey-Wilson g-derivative instead of the ordinary derivative. We have not explicitly 
worked out the limit transition g f 1 of the results, but by the set-up it is clear that the formal limit 
gives back many of the results of EH], EH- 

The contents of the paper are as follows. In Section E] we fix notation regarding matrix-valued or¬ 
thogonal polynomials. In Section [3] the notation for quantised universal enveloping algebras is recalled. 
Section 0] states all the main results of this paper. It introduces the quantum symmetric pair explicitly. 
Using the representations of the quantised universal enveloping algebra and the coideal subalgebra, the 
matrix-valued polynomials are introduced. We continue to give explicit information on the orthogonal¬ 
ity relations, three-term recurrence relations, g-difference operators, the commutant of the weight, the 
LDU-decomposition of the weight, the decoupling of the g-difference equations and the link to scalar¬ 
valued orthogonal polynomials from the g-Askey scheme. The proofs of the statements of Section |T] 
occupy the rest of the paper. In Section [S] the main properties derivable from the quantum group set-up 
are derived, and we discuss in Appendix |A] the precise relation of the branching rule for this quantum 
symmetric pair and the standard Clebsch-Gordan decomposition. In Section [5] we continue the study 
of the orthogonality relations, in which we make the weight explicit. This requires several identities 
involving basic hypergeometric series, whose proofs we relegate to Appendix [BJ Section [7] studies the 
consequences of the explicit form of the matrix-valued g-difference operators of Askey-Wilson type to 
which the matrix-valued orthogonal polynomials are eigenfunctions. 

In preparing this paper we have used computer algebra in order to verify the statements up to certain 
size of the matrix and up to certain degree of the polynomial in order to eliminate errors and typos. 
Note however, that all proofs are direct and do not use computer algebra. A computer algebra package 
used for this purpose can be found on the homepage of the second authorQ 

The convention on notation follows Kolb m for quantised universal enveloping algebras and right 
coideal subalgebras and we follow Gasper and Rahman [16] for the convention on basic hypergeometric 
series and we assume 0 < g < I. 
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2. Matrix-VALUED orthogonal polynomials 


In this section we fix notation and give a short background to matrix-valued orthogonal polynomials, 
which were originally introduced by Krein in the forties, see e.g. references in General references 

for this section are [5], [10], [18], and references given there. 

Assume that we have a matrix-valued function W: [a, &] —>■ M 2 e+i{C), 2£ -|- 1 S N, a < 6, so that 
W(x) >0 for X G [a, b] almost everywhere. We use the notation A > 0 to denote a strictly positive 
definite matrix. Moreover, we assume that all moments exist, where integration of a matrix-valued 
function means that each matrix entry is separately integrated. In particular, the integrals are matrices 
in M 2 f+i(C). It then follows that for matrix-valued polynomials P,Q G M 2 e+i{C)[x] the integral 

(2.1) {P,Q) = f P{x)*W{x)Q{x)dxGM 2 e+i{C). 

J a 

exists. This gives a matrix-valued inner product on the space M 2 i+i{C){x\ of matrix-valued polynomials, 
satisfying 

(P, Q) = (Q, P)\ (P, QA + RB) = (P, Q)A + (P, R)B, 

(P,P) = 0 G M 2 e+i (C) ^ Pix)=0G M 2 t+i (C) VX 

for all P,Q,R G M 2 e+i{C)[x] and A,B G M 2 e+i{C). More general matrix-valued measures can be 
considered cn], but for this paper the above set-up suffices. 

A matrix-valued polynomial P(x) = X(r=o i P^ ^ M 2 ^+i(C) is of degree n if the leading coefficient 
P" is non-zero. Given a weight W, there exists a family of matrix-valued polynomials (Pn)neN so that 
Pn is a matrix-valued polynomial of degree n and 

j-b 

( 2 . 2 ) / {Pn{x))*W{x)Pm{x)dx = 5n,mGn, 

J a 

where Gn > 0. Moreover, the leading coefficient of P„ is non-singular. Any other family of polynomials 
{Qn) SO that Qn is a matrix-valued polynomial of degree n and {Qn,Qm) = 0 for n ^ m, satisfies 
Pn[x) = Qn{x)En for some non-singular P„ G M2t+i{€.) for all n G N. We call the matrix-valued 
polynomial Pn monic in case the leading coefficient is the identity matrix I. The polynomials Pn are 
called orthonormal in case the squared norm = I for all n G N in the orthogonality relations (1^ . 

The matrix-valued orthogonal polynomials P„ satisfy a matrix-valued three-term recurrence of the 
form 


(2.3) xP„(x) = P„+i(x)A„ -I- P„(x)P„ -I- P„_i(x)G„ 

for matrices G M 2 e+i{C) for all n G N. Note that in particular A„ is non-singular for all n. 

Gonversely, assuming P_i(x) = 0 (by convention) and fixing the constant polynomial Po(x) G M 2 e+i{C) 
we can generate the polynomials Pn from the recursion (12.3p . In case the polynomials are monic, 
the coefficient A„ = I for all n and Po(x) = I as the initial value. In general, the matrices satisfy 
Gn+iAn = GnBn = B*Gn, SO that in the monic case Cn = G“^jG„ for n > I. In case the 

polynomials are orthonormal, we have G„ = A *_2 and P„ Hermitian. 

Note that the matrix-valued ‘sesquilinear form ’(El]) is antilinear in the first entry of the inner product, 
which leads to a three-term recurrence of the form (12.31) where the multiplication by the constant matrices 
is from the right, see m for a discussion. 

In case a subspace V C is invariant for W{x) for all x, is also invariant for W{x) for 

all X. Let Lv'-V^ be the embedding of V into so that Py = lyty G M 2 £+i(C) is the 

corresponding orthogonal projection, then W{x)Py = PvW{x) for all x. Let P^ : [a, h\ —>• End(fo)[x] be 
the matrix-valued polynomial defined by P^ {x) = 6yP„(x)ty, where P„ are the monic matrix-valued 
orthogonal polynomials for the weight W. Then P)[ form a family of monic fo-endomorphism-valued 
orthogonal polynomials, and Pn{x) = P^ (x) © P^^(x). The same decomposition can be written down 
for the orthonormal polynomials. 

The projections on invariant subspaces are in the commutant ^-algebra {T G M 2 i+i{C) \ TW{x) = 
W{x)T Vx}. In case the commutant algebra is trivial, the matrix-valued orthogonal polynomials are 
irreducible. The primitive idempotents correspond to the minimal invariant subspaces, and hence they 
determine the decomposition of the matrix-valued orthogonal polynomials into irreducible cases. 

Note that we only consider orthogonal decompositions of the matrix-valued weights and the corre¬ 
sponding polynomials. Tirao m discusses more general non-orthogonal decompositions. 
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We denote by Ei^j S M 2 ^+i(C) the matrix with zeroes except at the (j,j)-the entry where it is 1. So 
for the corresponding standard basis {ek}kLo set EijCk = Sj^kSi- We usually use the basis {ek}kLo 
in describing the results for the matrix-valued orthogonal polynomials, but occasionally the basis is 
relabelled as is customary for the Z/fq(su(2))-representations of spin i. In the latter case we 

use superscripts to distinguish from the previous case; E^i,j, k G {—i, ■ ■ ■ ,£}■ 


3. Quantised universal enveloping algebra 


We recall the setting for quantised universal enveloping algebras and quantised function algebras, and 
this section is mainly meant to fix notation. The definitions can be found at various sources on quantum 
groups, such as the books m, m, m, and we follow Kolb Ell- 

Fix for the rest of this paper 0 < q < 1. The quantised universal enveloping algebra can be associated 
to any root datum, but we only need the simplest cases g = s[(2) and g = s[(2) © s[(2). The quantised 
universal enveloping algebra is the unital associative algebra generated by k, k~^, e, / subject to the 
relations 

(3.1) kk~^ = l = k~^k, ke = q^ek kf = q~'^fk, ef — fe=- — -—r, 

q-q 

where we follow the convention as in [571 §3]- For our purposes it is useful to extend the algebra with 
the roots of k and k~^, denoted by satisfying 

,1/2, -1/2 _ , _ , -1/2, 1/2 l1/2^1/2 _ u ^-l/2jL-l/2 _ , -1 

f\j fXi .X. F\j i\j • rXj tXj fXi ■ f\j rv r\j ■ 

(3.2) , , , 

^ ^ k^l^e = qek^l\ k^/^ f = q-^ fk^/^. 


The extended algebra is denoted by Uq{sl{2)), and it is a Hopf-algebra with comultiplication A, counit 
e, antipode S defined on the generators by 


A(e) = e©l + fc©e, A(/) = / © + 1 ©/, A(fc=^i/2) = fc=^i/2 ^ ;.±i/2^ 

S'(e) = -fc"^e, S{f) = -fk, 5'(fc=^^/^) = e(e) = 0 = e(/), e{k^^/‘^) = l. 


The Hopf-algebra has a *-structure defined on the generators by 

(fc±i/2)* = fc±i/2^ r=q-^k-^e. 

We denote the corresponding Hopf *-algebra by Wq(su(2)). Then the identification as Hopf *-algebras 
with [30], HH is given by (A, B,C, D) O 

The irreducible finite dimensional type 1 representations of the underlying *-algebra have been clas¬ 
sified. Here type 1 means that the spectrum of is contained in For each dimension 2£ + 1 of 
spin i G there is a representation in with orthonormal basis {el^, • • • , e^} and on 

which the action is given by 

(3.3) , 1 /- 

b^(p) = - q^-P+^)(q-^-P - q^+P), 

where : Uq{5u{2)) —End('H^) is the corresponding representation. Note that b^{p) = b^{l—p). Finally, 
recall that the centre Z{Uq{5u{2))) is generated by the Casimir element w, 
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(3.4) 


UJ = 


j ^k ^ + qk — 2 

(9-1 - q? 


/e = 


qk 1 + q ifc — 2 
(9-1 - q? 


e/, 



We use the notation ^^(g) to denote the Hopf ^-algebra Uq{su{2)) © Uq{su{2)) — Uq{su{2) © su(2)), 
where Ei, Fi, / = 1, 2, are the generators. For fixed i the relations (EU cLnd. liold Avith. 

(fc^/^, e,/) replaced by , Ei, Fi) and the generators with different index i commute. 

The tensor product of two Hopf ^-algebras is again a Hopf ^-algebra, where the maps on a simple tensor 
Xi © X 2 are given by, see e.g. [SI Ch. 4], 

A(Ai © A 2 ) = Ai3(Ai)A 24(A2), e(Ai © A 2 ) = £(Ai)£(A2), 

(3.5) 


^(Ai © A 2 ) = ^(Ai) © 5'(A2), (Ai © A 2 )* = X* © A, 


2 , 


where we use leg-numbering notation. 

The irreducible finite dimensional type 1 representations of Uq{Q) are labeled by (€i,Q) G x ^N, 
and the representations : llq{Q) -G End('H^i’^2), = qiE ^ qiE^ g^j-g obtained as the exterior 
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tensor product of the representations of spin and ^2 oiUq[su.{2)) . Here type 1 means that the spectrum 
of i = 1,2, is contained in 

We have used the notation A, e, S for the comultiplication, counit, and antipode for all Hopf algebras. 
From the context it should be clear which comultiplication, counit and antipode is meant. The corre¬ 
sponding dual Hopf *-algebra related to the quantised function algebra is not needed for the description 
of the results in Section |TJ and it will be recalled in Section oi 


4. Matrix-valued orthogonal polynomials related to the quantum analogue of 

(SU(2) X SU(2),diag) 


In this section we state the main results of the paper. First we introduce the specific quantum sym¬ 
metric pair, which is to be considered as the quantum analogue of a symmetric space GjK, in this case 
SU(2) X SU(2)/SU(2). Quantum symmetric spaces have been introduced and studied in detail by Letzter 
[31], [3S], [31], see also Kolb m, In particular, Letzter has shown that Macdonald polynomials occur as 
spherical functions on quantum symmetric pairs motivated by work by Koornwinder, Dijkhuizen, Noumi 
and others. In our case, B C Kqid) as in Definition l4.Il is the appropriate right coideal subalgebra. Using 
the explicit branching rules for of Theorem 14.31 we introduce matrix-valued spherical functions 

in Definition 14.41 To these matrix-valued spherical functions we associate matrix-valued polynomials in 
(jUj), and we spend the remainder of this section to describe properties of these matrix-valued polynomi¬ 
als. This includes the orthogonality relations, the three-term recurrence relation, and the matrix-valued 
polynomials as eigenfunctions of a commuting set of matrix-valued ^-difference equations of Askey-Wilson 
type. Moreover, we give two explicit descriptions of the matrix-valued weight function W, one in terms 
of spherical functions for this quantum symmetric pair, and one in terms of the LDU-decomposition. 
The LDU-decomposition gives the possibility to decouple the matrix-valued g-difference operator, and 
this leads to an explicit expression for the matrix entries of the matrix-valued orthogonal polynomials in 
terms of scalar-valued orthogonal polynomials from the g-Askey scheme in Theorem 14.171 

For the symmetric pair {G,K) = (SU(2) x SU(2), SU(2)), K = SU(2) corresponds to the fixed points 
of the Cartan involution 9 flipping the order of the pairs in G. The quantised universal enveloping 
algebra associated to G is Uq{Q) as introduced in Section |31 As the quantum analogue of K we take 
the right coideal subalgebra B C Uq{Q), i.e. B C Uq[Q) is an algebra satisfying A(H) d B ® Uq[Q), as 
in Definition oi Letzter [31 Section 7, (7.2)] has introduced the corresponding left coideal subalgebra, 
and we follow Kolb [27l §5] in using right coideal subalgebras for quantum symmetric pairs. Note that 
we have modified the generators slightly in order to have H* = H 2 . 


Definition 4.1. The right coideal subalgebra B C Gq{Q) is the subalgebra generated by ^ where 

K = KiK^^, and 

Hi = + qF^K-^'^K^^, + qF^R^^^ R-^^^. 


Remark 4.2. (i) H is a right coideal as follows from the general construction, see [23 Prop. 5.2]. It can 
be verified directly by checking it for the generators. Note = R^^G (g) R^'^G is immediate, 

and 

A(Hi) = Hi (g) {RxR^Y^G + if 1/2 ^ q-^{R^R^)-^GEi -f R-^G ^ qP^R-^G g h (g) W,(g) 

by a straightforward calculation. Since H 2 = B^, it also follows for H 2 , since R^^G are self-adjoint. The 
relations, cf. EZl Lemma 5.15], 

(4.1) RG^B^=qBiRG\ rG^ = q-^ B^R^^, [Hi,H 2 ] = 

q-q-^ 

hold in Uq{Q) as can also be checked directly. 

(ii) Let T: Uq{su{2)) — >• Uq{su{2)) be defined by 

4 -(fci/ 2 ) = fc- 1 / 2 ^ ^ ^ 1 / 2 ^ ^f(e) = g3/, ^^l(f)=q-^e, 

then 41 extends to an involutive *-algebra isomorphism. Consider the map 

i o ( 4 / (g) Id) o A: Uq{su{2)) Wq(0) 

where t is the algebra morphism mapping a; (g) y £ Uq{su{2)) (g) Uq{su{2)) to the corresponding element 
X 1 Y 2 £ hlq{Q) for X and y generators of Uq{su{2)). Then we see that YG i_>. qfYG i_>. Hi, and 

q~^k~^G^ 1 —^ H 2 under the map t o (41 g) Id) o A. In particular, the relations (14.11) follow. We conclude 
that B is isomorphic as a *-algebra to A(t7g(su(2)) C Uq{'Y) by the *-isomorphism 60 (41 g Id). 
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(iii) In particular, B = Uq{su{2)) as ^-algebras. So the irreducible type 1 representations of B are 
labeled by the spin (. G ^N. This can be made explicit hy : B ^ End(H^) and setting 

(4.2) = q~Pel, t\B2)el = b\p + l)e^+i 

with the notation of (|3.3I) . We use the same notation for these representations here and in (13.311 . since 
they correspond under the identification of B as Uq{su{2)). 

(iv) Let tr be the ^-algebra isomorphism on Uq{Q) — Uq{5u{2)) ®Uq{5u.{2)) by flipping the order in the 

tensor product, or equivalently by flipping the subscripts 10 2. Then cr: O is an involution Bi o 
B 2 j K o K~^. On the level of representations of Wq(g) and B it follows {X)P, 

X G Uq{g), where P: O jg the flip, and t\aiY)) = {J^)*t^{Y)j\ Y G B, where ^ 

n^, J^: e^O eip. 


Theorem 4.3. The finite dimensional representation ofUq{g) restricted to B decomposes multi- 

plieity-free into irreducible representations t^ of B; 






0 




0 




Moreover, with respect to the orthonormal bases {ep}p^_^ of TL^ and 0 for 

the B-intertwiner TL^ ^ is given by 




E 

i=-ti 3 = -l 


i2 

E 


cf 


e ei 
e/ 




where are Clebsch-Gordan coefficients satisfying = 0 ifi — jf^p. 

The proof of Theorem 14.31 is a reduction to the well-known Clebsch-Gordan decomposition for the 
quantised universal enveloping algebra Uq{su{2)), see e.g. [5], [H], using Remark 14.21 The proof is 
presented in Appendix 1X1 In particular, {j3g^ Pti A identity on . Note that 




A,A 


-H" 


0 e 


A 

"2 




'-'ni,n2,p 

p——£ 


e 


l 

p- 


In general the decomposition of an irreducible representation restricted to a right coideal subalgebra 
seems a difficult problem. In this particular case we can reduce to the Clebsch-Gordan decomposition, 
and yet another known special case is by Oblomkov and Stokman [351 Prop. 1.15], but in general this is 
an open problem. 

In particular, for fixed £ G ^N, we have [t^^’^^ls- t^] = 1 if and only if 

(4.3) (£i,£ 2) e ^Nx ^N, \£i-£2\<£<£i+£2, £i + £2-£€^- 

We use the reparametrisation of (14.3p by 

(4.4) C = e':Nx{0,l,...,2£}^iNxiN, ^(n, fc) = ( ^ 



see also Figure [Hand [23l Fig. 1, 2]. In case £ = 0, we have = e|e, where e is the counit of Uq{Q), is 
the trivial representation of B, and the condition (ESI) gives £i = £2 and f^{n,0) = {^n, ^n). 

With these preparations we can introduce the matrix-valued spherical functions associated to a fixed 
representation t^ of B, where we use the notation of Theorem 14.31 


Definition 4.4. Fix £ G and assume (£i,£ 2 ) G x so that t^] = 1. The spherical 

function of type £ associated to (£i,£ 2 ) is defined by 

■■ ^ End(?£0, ^ ^ o o 

Remark 4.5. (i) Note that the requirement on (£i,£ 2 ) in Definition 14.41 corresponds to the condition 
(031). Since is a S-intertwiner, we have 

(4.5) <,,jAzr) = t^(AX,,jz)£^(r), yx,YGB,'iZGUqiQ). 
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(ii) Note that the condition (14.31) is symmetric in l\ and li. With the notation of Remark 14.21' ivi 
we have 4)^^ ^ £ ^ 9 ( 0 )- This follows from which is a 

consequence of (TOI . 

In case £ = 0, = C, we need £1 = £ 2 - Then are linear maps Uq{Q) —>■ C. In particular, $0 g 

equals the counit e, and the spherical function ip = \{q~^ + q )^\/2 1/2 scalar-valued linear map on 
Uq{Q). The elements n /2 written as a multiple of Un{p>), where Un denotes the Chebyshev 

polynomial of the second kind of degree n, see Proposition 15.31 This statement can be considered as 
a special case of Theorems 14.81 but we need the identification with the Chebyshev polynomials in the 
spherical case £ = 0 in order to obtain the weight function in Theorem 14.81 Proposition 15.31 will follow 
from Theorem IMl The identification of the spherical functions for £ = 0 with Chebyshev polynomials 
corresponds to the classical case, since the spherical functions on G x G/G are the characters on G and 
the characters on SU(2) are Chebyshev polynomials of the second kind, as the simplest case of the Weyl 
character formula. It also corresponds to the computation of the characters on the quantum SU(2) group 
by Woronowicz 07], since the characters are identified with Chebyshev polynomials as well. 

Next Theorem l4.6l gives the possibility to associate polynomials in (p to spherical functions of Definition 
14.41 Theorem 14.61 essentially follows from the tensor product decomposition of representations of Gg(g), 
which in turn follows from tensor product decomposition for Uq(su(2)), and some explicit knowledge of 
Clebsch-Gordan coefficients. 

Theorem 4.6. Fix £ £ and let (£i,£ 2 ) £ x satisfy (14.31) . then for constants Aij we have 

= X! ^ 1 / 2 , 1/2 0. 

i,3 = ±l/2 

In order to interpret the result of Theorem l4.6l we evaluate both sides at an arbitrary X £ Gg(g). The 
right hand side is a linear combination of linear maps from to itself after evaluating at X. For the 
left hand side we use the pairing of Hopf algebras, so that multiplication and comultiplication are dual 
to each other and the left hand side has to be interpreted as 

(4.6) a) (X) = Y. <.^ 2 (^( 2 )) e End(H^), 

(X) 

which is a linear combination of linear maps from T-L^ to itself, using A(X) = X)(x) "^( 1 ) ® "^(2)- The 
convention in Theorem 14.61 is that Ai^j is zero in case (£1 -|- ^,£2 -I-j) does not satisfy (14.311 . The proof of 
Theorem 14.61 can be found in Section O 

Since B is right coideal subalgebra we see that the left hand side of Theorem 14.61 has the same 
transformation behaviour as (1131) . Indeed, for A £ and Y £ Uq{Q) we have 

(p^Y^){XY)= Y ‘£’(^(i)T(i)X.^,(^(2)T(2)) 

{x),(X) 

(4 7) ^ X! ^(^(i))y’(T(i)) .^(-’^(2)^(2)) = y^y?(F(i)) 4>^,,f,(y^£(A(i))A(2)y(2)) 

{X),(Y) (Y) (X) 

= Yf{Y(i)) <,,,(AY(2)) =YF{Y(i))t\X)^Y2(Y(2)) = t\X){p<^Y^YY), 

{Y) {Y) 

where we have used that £ S by the right coideal property and (14.51) for p = 4>°y2 1/2 in f^e second 
equality, and the counit axiom e(A(i))A( 2 ) = X in the fourth equality and then (14.51) for 

and the fact that p{Y(^i^) is a scalar. Similarly, the invariance property from the right can be proved. 
Theorem 14.61 leads to polynomials in p by iterating the result and using that Ai/ 2 , 1/2 is non-zero. 

Corollary 4.7. There exist 2£ + 1 polynomials 0 < k < 2£, of degree at most n so that 
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= XI ^n,m(V5) 4>4(o,fc), U £ N, 0 < m < 2£. 

k=0 

The aim of the paper is to show that the polynomials give rise to matrix-valued orthogonal 
polynomials. Put 

P„ = P^£End(H0N {Pn\j=^p 


(4.8) 


0<z,j<2£, 




Figure 1. The spherical functions when ^ = 2 and interpretation of 5/2 

in terms of the matrix-valued spherical functions. The reparametrization ^ is depicted. 


where the matrix-valued polynomials are taken with respect to the relabeled standard basis Cp = ep_^, 
p S {0,1, • • • , 2£}, so that P„ = j=o ® ■ From Corollary 14.121 or Theorem 14.171 we see that 

the polynomial has real coefficients. The case £ = 0 corresponds to a three-term recurrence rela¬ 
tion for (scalar-valued) orthogonal polynomials, and then the polynomials coincide with the Chebyshev 
polynomials [/„ viewed as a subclass of Askey-Wilson polynomials [U (2.18)], see Proposition 15.31 

We show that the matrix-valued polynomials {Pn)'^=o are orthogonal with respect to an explicit 
matrix-valued weight function W, see Theorem 14.81 arising from the Schur orthogonality relations. The 
expansion of the entries of weight function in terms of Chebyshev polynomials is given by quantum group 
theoretic considerations except for the calculation of the coefficients in this expansion. The matrix-valued 
orthogonal polynomials satisfy a matrix-valued three-term recurrence relation as follows from Theorem 
14.61 which in turn is a consequence of the decomposition of tensor product representations of lAq{Q). 
However, in order to determine the matrix coefficients in the matrix-valued three-term recurrence we use 
analytic methods. The existence of two Casimir elements in Uq{Q) leads to the matrix-valued orthogonal 
polynomials being eigenfunctions of two commuting matrix-valued ^-difference operators, see m for the 
group case. This extends Letzter |35j to the matrix-valued set-up for this particular case. The g-difFerence 
operators are the key to determining the entries of the matrix-valued orthogonal polynomials explicitly 
in terms of scalar-valued orthogonal polynomials from the g-Askey-scheme [26) . namely the continuous 
g-ultraspherical polynomials and the g-Racah polynomials. In this deduction the LDU-decomposition of 
the matrix-valued weight function W is essential, since the conjugation with L allows us to decouple the 
matrix-valued g-difference operator. 

In the remainder of Section |4] we state these results explicitly, and we present the proofs in the 
remaining sections. First we give the main statements which essentially follow from the quantum group 
theoretic set-up, except for explicit calculations, and these are Theorems 14.81 ICTl 14.131 The remaining 
Theorems 14.15114.171 are obtained using scalar orthogonal polynomials from the g-analogue of the Askey 
scheme m and transformation and summation formulas for basic hypergeometric series m- 

We start by stating that the matrix-valued polynomials {Pn)^=o introduced in (|4.8I1 are orthogonal 
with the conventions of Section [2] The orthogonality relations of Theorem 14.81 are due to the Schur 
orthogonality relations. The expansion of the entries of the weight function in terms of Chebyshev 
polynomials follows from the fact that the entries are spherical functions, i.e. correspond to the case 
= 0 so that they are polynomial in ip. The non-zero entries follow by considering tensor product 
decompositions, but the explicit values for the coefficients at(rn,n) in Theorem 14.81 require summation 
and transformation formulae for basic hypergeometric series. 

Theorem 4.8. The polynomials (Prt))]ffio defined in (14. 8|) form a family of matrix-valued orthogonal 
polynomials, so that Pn is of degree n with non-singular leading coefficient. The orthogonality for the 
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matrix-valued polynomials {Pn)n>o is given by 

-[ Pm{x)*W{x)Pn{x)\/l - X'^dx = Gm5m,n: 
TT 7-1 

where the squared norm matrix Gm is diagonal; 

(1 5«+2)2 


{Gn)iJ — 


.2n-2e 


_ q2n+2i+2'^(J^ _ qU-2i+2n+2y 

Moreover, for 0 < m < n < 2£ the weight matrix is given by 

n 

^ ^ OLtijn, Ti) U'm-\-n—2t(,x\ 
t=0 

where 

n, (m n') — „2ri.(2^+l)-n^-(4£+3)t+t=^-2^+m 1 - \ q^)2i-n{q^ \ q^)r 

l-g“ {q^-,q^)2i 

(_l)m-t(^2m-«. (^4^+4-2t. q2y 


(^2m+4.g2)^_^ 




and W{x)m,n = W{x)n,m if m > u. 


The proof of Theorem 14.81 proceeds in steps. First we study explicitly the case i = 0, motivated by 
works by Koornwinder [30], Letzter [34], [35], [36], and others. Secondly, we show that taking traces of a 
matrix-valued spherical function of type i associated to times the adjoint of a spherical function 

of type i associated to {£[,£' 2 ) gives, up to an action by an invertible group-like element of Uq{Q), a 
polynomial in the generator for the case £ = 0. Then the explicit expression of the Haar functional on 
this polynomial algebra, stated in Lemma |5.41 gives the matrix-valued orthogonality relations. Finally, 
the explicit expression for the weight is obtained by analysing the explicit expression of W in terms of 
the matrix entries of the intertwiners in case €1 4-^2 = £■ These matrix entries are Clebsch-Gordan 

coefficients. 

The leading coefficient of P„ can be calculated explicitly from the proof of Theorem 14.81 


Corollary 4.9. 


The leading coeffieient of Pn is a non-singular diagonal matrix; 


(„2i+2 M-2i+2 2\ 

Ay I'i )r 


(g2,g«+4.^2)^ 


The weight W is not irreducible, see Section [2] but splits into two irreducible block matrices. The 
symmetry J of the weight function of Theorem l4.8l is essentially a consequence of Remark ld.sr iib but we 
need the explicit expression of the weight in order to prove that the commutant algebra is not bigger. 


Proposition 4.10. The commutant algebra 

{W{x) I X G [-1,1]}' = {r G EndCH^) I W{x)Y = YW{x),\/x G (-1,1)}, 

is spanned by I and J, where J: Cp e 2 i-p, p G {0, • • • , 2£}, is a self-adjoint involution. Then 
JPn{x)J = Pn{x) and JGn(x)J = Gn{x). Moreover, the weight W decomposes into two irreducible block 
matrices W+ and W-, where IV+, respectively W-, acts in the +l-eigenspace, respectively —1-eigenspace, 
of J. So for P+ + P_ = I, where P+, P- are the orthogonal self-adjoint projections P+ = \{I + J), 
P- = — we have that W+, respectively W-, corresponds to P+W{x)P+, respectively P-W{x)P-, 

restricted to the +l-eigenspace, respectively —1-eigenspace, of J. 

The special cases ior £ = ^ and £ = 1 are given at the end of this section. In particular, we identify 
all scalar-valued orthogonal polynomials occurring in this framework explicitly in terms of Askey-Wilson 
polynomials. 

Theorem 14.61 can be used to find a three-term recurrence relation for the matrix-valued orthogonal 
polynomials P„, cf. Section [S] so the underlying tensor product decompositions provide the three-term 
recurrence relation. However, the resulting expressions for the entries of the coefficients of the matrices 
are rather complicated expressions in terms of Clebsch-Gordan coefficients. For the corresponding matrix¬ 
valued monic polynomials Qn{x) = P„(x)lc(P„)“^, see Corollary 14.91 for the explicit expression for the 
leading coefficient, we can derive a simple expression for the matrices in the three-term recurrence relation 
once we have obtained more explicit expressions for the matrix entries of Qn- This is obtained in Section 
[7] using an explicit link of the matrix entries to scalar orthogonal polynomials in the g-Askey scheme. 
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Theorem 4.11. The monic matrix-valued orthogonal polynomials {Qn)n>Q satisfy the three-term recur¬ 
rence relation 

where Q-i(x) = 0, Qo(x) = I and 

1 g2n+l — g4^+2"+2)2 

( 2 (1 _ q2i+2n^(^^ _ g4£+2n-2i^^2 — g2r!,+2i+2^^]^ _ g4£-2i+2n+2^ 

^q2n+l _ ^2n)2('^ _ ^«+2n+2)2 

" 2^ 2 (1 — q2n-\-2i^(^l _ qie-^-2n-2i'^(^l _ g2n-t-2i-t-2^ _ ^«-2*+2n+2 ^’ 

1 (l-g2")2(i_^4£+2„+2)2 

4 (1 _ q2n-\-2i'^(^l _ q4e-\-2n-2i'^(^l _ ^2n+2i+2)('2^ _ ^4^-2j+2n+2)'^*4' 

Note that —>• 0, as n —>■ oo. 

The three-term recurrence relation for the matrix-valued orthogonal polynomials P„ is given in Corol- 
larv l4.12l which follows from Theorem ll.lll since G„+iA„ = lc(P„+i)*lc(P„), GnBn = lc(P„)*X„lc(P„), 
Gn-iCn = lc(P„_i)*F„lc(P„). For future reference we give the explicit expressions in Corollary 14. 121 

Corollary 4.12. The matrix-valued orthogonal polynomials {Pn)n>o satisfy the three term recurrence 
relation 


11 - 




2=C 






z=0 


xPn(^^^ — -fn+l(^)-^n “1“ -fn(^)-^n H” -fn— 


where P-i{x) = 0, Pq{x) = I and 


, ^ 1 (l-g2"-H2)n_ 4^-S2„-h4) 

An = Y. — .. ... 


2=0 


2q (1 - g2z+2n+2)^X - ^4^-22+271+2)^^^^^’ 


_2^+2n+l (i_g4f-2.)(i_^2.+2) 

" 2 (1 — g4£-l-2n-2j)('X _ g2n-|-2i-|-4^ *>*+1 


le 


, _ V- /V- ^ _ / p 

2 (1 - g2n-|-2i^('2^ _ ^4£-2i-|-2n-|-4) *>*-!> 

(l_g2»)(i_^4mn+2) 

^ (y /-I 9^_L9-)_l_0\ /1 /I ^-LO-ki — Qo-LO \ ' 


2n+l 


(1-q2*)(1_^4£-22+2) 


2=0 


2 _ ^271+22+2)^^ _ g4£+272—22+2) ’ 


Note that the case i = 0 gives a three-term recurrence relation that can be solved in terms of the 
Chebyshev polynomials, see Proposition 15.31 

In the group case, the spherical functions are eigenfunctions of it'-invariant differential operators on 
G/K, see e.g. [5], [IS]- For matrix-valued spherical functions this is also the case, see mi, and this has 
been exploited in the special cases studied in m, m, m- In the quantum group case the action of the 
Casimir operator gives rise to a g-difference operator for the corresponding spherical functions, see [35] , 
The first occurrence of an Askey-Wilson g-difference operator, see a, m, m. in this context is due 
to Koornwinder [30| . For the matrix-valued orthogonal polynomials we have a matrix-valued analogue 
of the Askey-Wilson g-difference operator, as given in Theorem 14.131 We obtain two of these operators, 
one arising from the Casimir operator for Uq{su{2)) in the first leg of Gg(g) and one from the Uq{su{2)) 
Casimir operator of the second leg. This is related to a kind of Cartan decomposition of l{q{g), cf. (14.5L 
which, however, does not exist in general for quantised universal enveloping algebras. We can still resolve 
this problem using techniques based on [5J §2], see the first part of the proof in Section [SJ The proof of 
Theorem 14.131 is completed in Section [T] 


Theorem 4.13. Define two matrix-valued q-difference operators by 

Di — A^iii^Z^rjq -\- — ^ 7 I — 1 - 5 ^; 


where the multiplication by the matrix-valued functions M.i{z) and A4i(z is from the left, and where 
Pq is the shift operator defined by {pqf){z) = f{qz), f{z) = f{p,{z)), where x = pL{z) = \{z + z~^). The 
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matrix-valued function Aii is given by 


2i-l 


Mi[z)=- 


*(1 —z 


2i 


2=0 




- + n- 


® (1 — g^*+^z^) 


i=0 


(1 _ ^2)2 (1 _ ^2) 


-^2,2 5 


and Ai 2 {z) = JA4i{z)J, where Jcp = e^n-p- The matrix-valued orthogonal polynomials Pn are eigen¬ 
functions for the operators Di with eigenvalue matrices given by A„(i) such that DiPn = Pn^nii) o-nd 
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„-3-n-l I „j+n+l 

A»(1) = E .-1 ^2 -A„(2) = JA„(1)J. 


j=o 


(g-i - g)2 


Explicitly 


{DiPn){fJ.{z)) = Mi{z)(pqP„)(z) + Mi(z ^ P„) (z) = Pn{p{z))An{i), 


where pq and Pq-i are applied entry-wise to the matrix-valued orthogonal polynomials Pn- 

Theorem 14.131 shows that JDiJ = D 2 , since J is constant. In particular, Di + D 2 commutes with J 
and reduces to a g-difference operator for the matrix-valued orthogonal polynomials associated with the 
weight W+ or W -, see Proposition 14.101 Similarly, Di — D 2 anticommutes with J. 

Note that for any matrix-valued polynomial P, the expression Mi{z)P{qz) -\- Mi{z~^)P{z/q) is sym¬ 
metric in z o z“^, and hence again a function in a: = fJ,{z). The case £ = 0 corresponds to only one 
g-difference operator, which we rewrite as 

(4.9) P" = (<?■” + 

The relation (14.9p has the Chebyshev polynomials Un{x) = (g”'*'^; g)^^Pn(a;; g, —g, g^^^, —g^^^lg), rewrit¬ 
ten as Askey-Wilson polynomials [H (2.18)], as solutions, see [551 §14-1], [IB §7.7], [551 Ch. 15-16]. In 
particular, we consider the operators of Theorem 14.131 as matrix-valued analogues of the Askey-Wilson 
operator, see Askey and Wilson a , or a, m. m- 


Corollary 4.14. The q-difference operators Di and D 2 are symmetric with respect to the matrix-valued 
weight W, i.e. for all matrix-valued polynomials P, Q, we have 



{{D^P){x))*W{x)Q{x) dx 



{P{x))*W{x){D^Q){x) dx, 


1 = 1,2. 


By [IHl §2] it suffices to check Corollary 14.141 for P = P„, Q = Pm, so that by Theorem 14.131 and 
Theorem 14.81 we need to check that An{i)*GnSm,n = GnAm{i)Sm,n, which is true since the matrices 
involved are real and diagonal. 

In order to study the matrix-valued orthogonal polynomials and the weight function in more detail 
we need the continuous g-ultraspherical polynomials [5] Chapter 2], [T5], [IB Chapter 20], [55] : 


(4.10) 


Cn{x;l3\q) 


n 


E 


(/3;g)r(/3;g)«-r ^i(n-2r)0 
{q;q)r{q;q)n-r 


X = COS 9. 


The continuous g-ultraspherical polynomials are orthogonal polynomials for |,5| < 1. The orthogonality 
measure is a positive measure in case 0 < g < 1 and /3 real with |/3| < 1. Explicitly 


(4.11) 


C'„(a:;/3|g)C'm(a:;/3|g) 


w{cos9; fi\q) = 


dx — dnm, 57r 


w{x;l3 I g) 
y/1 — x"^ 
(e^^^,e~^^^;g)oo 

(,5e2*^,/3e-2*®;g)oo' 


(/3,/3g;g)oo (/3^;g)n 1-/3 
(/3^, 9 ; 9)00 (g;g)n l-^g" 


Note that in the special case fd = q^'^^, k gN, the weight function is polynomial in a; = cos0, and 
(4.12) ri;(cos 6»; g^+'= |g) = 4(1 - cos^ 9) (ge^*^ q)k 


We use the continuous g-ultraspherical polynomials (14.101) for any /3 G C. In particular, for fd = q~^ 
with fc £ N the sum in (]4.10p is restricted to n — k < r < k, and in particular Gn{x', q~^;q) = 0 in case 
n — k > k. With this convention we can now describe the LDU-decomposition of the weight matrix, and 
state the inverse of the unipotent lower triangular matrix L in Theorem 14.151 
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Theorem 4.15. The matrix-valued weight W as in Theorem \4-8\ has the following LDU-decomposition: 

W{x) = L{x)T{x)L{xY , xS[—1,1], 

where L: [—1,1] —>■ M 2 ^+i(C) is the unipotent lower triangular matrix 

m-k 


L{x)mk = q 




{q'^]q^)m+k+i{q^-,q^)k 

and T: [—1,1] —> M 2 f+i(C) is the diagonal matrix, 0 < k < 2£, 


0 < k < m < 2i, 


Tix)kk = Ckii) = 


1-0^2 ’ 4 ( 9 ^; 92 ) 2 ,^i(g 2 ;g 2)2 

The inverse of L is given by 


2£. 


(g > g )y- ct-„(x; ,-^V). 0<n<L 

(g^;r)2fe(r;r)n 


Note that T, L and L~^ are matrix-valued polynomials, which is clear from the explicit expression 
and (14.121) . It is remarkable that the LDU-decomposition is for arbitrary size 2£-f 1, but that there is no 
dependence of L on the spin £ and that the dependence of T on the spin £ is only in the constants Ck {£) ■ 
We prove the first part of Theorem l4.15l in Section [S] The proof of Theorem l4.15l is analytic in nature, 
and a quantum group theoretic proof would be desirable. The statement on the inverse of L{x) is taken 
from [T], where the inverse of a lower triangular matrix with matrix entries continuous q-ultraspherical 
polynomials in a more general situation is derived. The inverse L~^ is derived in [TJ Example 4.2]. The 
inverse of L in the limit case g t 1 '''^as derived by Cagliero and Koornwinder [7] , and the proof of [T] is 
of a different nature than the proof presented in [7]. 

Theorem l4.15l shows that det(tT (x)) is the product of the diagonal entries of r(a:). Since all coefficients 
Ck{£) > 0 and the weight functions are positive, we obtain Corollary 14. 161 


Corollary 4.16. The matrix-valued weight W(x) is strictly positive definite for x G [—1,1]. In particular, 
the matrix-valued weight W{x)'/l — x"^ of Theorem\4-8\ is strictly positive definite for x € (—1,1). 


Using the lower triangular matrix L of the LDU-decomposition of Theorem l4.15l we are able te decouple 
Di of Theorem 14.131 after conjugation with L*{x). We get a scalar g-difference equation for each of the 
matrix entries of V{x)Pn{x), which is solved by continuous g-ultraspherical polynomials up to a constant. 
Since we have yet another matrix-valued q-difference operator for L*(x)Pn{x), namely nD 2 {V)~^ with 
D 2 as in Theorem 14.131 we get a relation for the constants involved. This relation turns out to be a 
three-term recurrence relation along columns, which can be identified with the three-term recurrence for 
g-Racah polynomials. Finally, use to obtain an explicit expression for the matrix entries of 

the matrix-valued orthogonal polynomials of Theorem 14.171 

Before stating Theorem 14.171 recall that the g-Racah polynomials, see e.g. [161 §7-2], |20l §15.6], [26l 
§14.2], are defined by 


(4.13) 


R„(^(x);a,/3,7,5;g) = 4‘/?3 


g-",a/3g”+i,g-"^,7(5g"^+i 

ag,/ 3 ^g, 7 <? 



nG{0,l,2,...,7V}, 


where iV G N, ix{x) = q ^ -\- and so that one of the conditions aq = q ^, or fidq = q ^, or 

7 g = q~^ holds. 


Theorem 4.17. For 0 < i,j < 2£ we have 
2t 

Pn{x)ij = y^(-l)'=^g"+(2fc+l)(fc-b+t(2fc-l-l)-l-2fc(2^+»+l) 


-k^ 


k=i 


(g^g^)fc(g^;g^)fc-^^ (g ^^g (g^;g^)n-^J--fc 

{q^;q^)2k{q^;q^fi {qfi q^^+^; q^)k {q^'^+^-, q^)n+j-k 

X Rki^l{J);l,l,q-^^-^^-^,q-'^‘^-^;q'^)Ck-^{x■,q-^'^\q‘^)Cn+J-kix■,q^^+‘^\q'^). 

Note that the left hand side is a polynomial of degree at most n, whereas the right hand side is of 
degree n -\- j — i. In particular, for j > i the leading coefficient of the right hand side of Theorem 14.171 
has to vanish, leading to Corollary 14.181 
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Corollary 4.18. With the notation of Theorem\4.17\ we have for j > i 


2^ / 2 2\ ( 2 2\ / —4£ 

E /_-| \fc (2fe+l)(fe-i)+j(2fc+l)+2fc(2f+n+l)-fc^ W i g )k\<l ] <1 )k+i jg _ ,JI_ 

^ {q^;q^)2k - 


-U „-2j-2n.g2)^ 


k—i 


(<72,g4^+4;g2). 


{q‘^^+^\q^)n+3-k 


Rk{t^{j)-A,l,q 


— 2n—2j—2 


-4i-2 


;g") = o. 


By evaluating Corollary 14.71 at 1 G htqid) we obtain Corollary 14.191 which is not clear from Theorem 
Corollary 4.19. For all m G {0, ■ ■ ■ , 2g} we have Yll^=o{Rri{\{q + q~^)))k,m = 1- 


4.1. Examples. We end this section by specialising the results for low-dimensional cases. The case 
£ = 0 reduces to the Chebyshev polynomials Un{x) of the second kind as observed following Theorem 
14.131 This is proved in Proposition 15.31 which is required for the proofs of the general statements of 
Section 01 


4.1.1. Example: £ = ^. In this case we work with 2x2 matrices. By Proposition 14.10| we know that 
the weight is block-diagonal, so that in this case we have an orthogonal decomposition to scalar-valued 
orthogonal polynomials. To be explicit, the matrix-valued weight W is given by 




0 

-2x + {q + q~^) 



X G [-1,1]. 


In this case, see Section 01 the polynomials diagonalise since the leading coefficient is diagonalised by 
conjugation with the orthogonal matrix Y, and we write YPn{x)Y* = diag(p+(a;),p“(x)). Then we can 
identify by any of the results given in this section, and we do this using the three-term recurrence 
relation of Corollary 14.121 After conjugation, the three-term recurrence relation for is given by 


XPn {x) = 


1 (l-g2"+6) 

2<7 (1 - + 


2n+l 


(i-a 


, 2\2 


2 (l-Q2n+2)(l_^2n+4) 


Pn (x) 


q jl-qn 

2 (1 - g2"-h2) 


Pt-iix), 


and the three-term recurrence relation for p^ is obtained by substituting x i—>■ —x into the three-term 
recurrence relation for p+. The explicit expressions for and pf are given in terms of continuous 

g-Jacobi polynomials P^'^\x\q) for (a,/3) = (i, |), see [H §4], 021 §14.10]. Note that xjg^) = 

(—l)^q~'^Pn‘^'^\x\q^), see [161 Exercise 7.32.(ii)]. So we obtain 


Pn (x) = 


(I_g2n-H4) ('q2n+6.^2)^ 


- P 

J 1 


(1 1 ) 


ix\q ), 


Pn (X) = 


(l-g4) (q^,-q3,-q4;q^), 

(I_g2n-h4) (^2n-H6.q2)^ 


P 


(1 1 ) 
2 ’ 2 '' 




which is a g-analogue of [231 §8.2]. Moreover, writing down the conjugation of the g-difference operator 
Di + D 2 of Theorem |4T3] for the case £ = ^ for the conjugated polynomials gives back the Askey-Wilson 

g-difference for the continuous g-Jacobi polynomials Pn°‘’^\x\q) for (a,/3) = (^, |) and (|, ^). Working 
out the eigenvalue equation for Di — D 2 gives a simple g-analogue of the contiguous relations of 031 
p. 5708]. 


4.1.2. Example: £ = 1. For £ = 1 we work with 3x3 matrices. By Proposition 14.101 we can block- 
diagonalise the matrix-valued weight: 


YW{x)Y* 



0 ^ 

W-ix)J ’ 


0 P 

Y = ^V2 I 0 72 0 

\-l 0 1; 


X G [-1,1], 


where W+ is a 2 x 2 matrix-valued weight and W- is a scalar-valued weight. Explicitly, 

.2 , „4^ 


W+{x) = 


4a;2-h(g2-hg 7 272g ^ 

(1 + r) 

4g2 


\ (1 + r) 


(l + g2)"‘ 


+ (7 + 7-7 


W-{x) = -47 + {q^+ 2 + g- 7 . 
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The polynomials diagonalise by YPn{x)Y* = diag(P+(x),p“(x)), where is a 2 x 2 matrix-valued 
polynomial and p~ is a scalar-valued polynomial. Conjugating Corollary 14.121 the three-term recurrence 
relations for and p~ are 


where 


xPi{x) = P++i(a:)A„ -h P+(a;)P„ -b P^_^{x)Cn, 

, 1 (l-g^»+8) _ , , , g (1-g^") - 

XPn [x) _ g2n+6'jPn+l(x) 2 (1 — q2n+2^Pn-l 


(x), 



/(I 

_ g2n+8'^ 

0 ^ 
(l-g2"+8)(l-g2"+2) 


/ (l-g2") 

0 ^ 
(1 - g2")(l - g2"+6) 

2g 

(1 

_ q2n+6^ 

0 

II 

C) 

(1 - g2"+2) 

1 " 


V 

(1 — q2n-|-4^2 J 


(1 - g2"+4)2 y 


Br 




^ 2 n+l 


0 

(l-g2)(i-54) 


2 n-n (l-g^)(l-g'^) \ 


(1 — q2n-|-4^2 




The scalar-valued polynomial can be identified with the continuous g-ultraspherical polynomials; 

r)-(x] — g" (1 ~ g )(1 ~ g ) ^ (x- 

Pn [ X ) - q _ ^2n-|-2)('X - g2n-|-6) '-'^\ X , q \q ). 

The 2x2 matrix-valued polynomials are solutions to the matrix-valued g-difference equation DP^^ = 
P+A„. Here D — M{z)rig+M{z~^)qq-i is the restriction of the conjugated D 1 +D 2 to the -bl-eigenspace 
of J. The explicit expressions for M{z) and A„ are 


M{z) = 


T-l 


(l-g 2 )(l-z 2 ) 


(l-g 2 )^ 


V 


-V2g(l 


q^z^) 

q^)z 


2 g 


—\/ 2 g^z 

(1 - q^z^) 

( 1 -g^) y 


A 


n 


/ _„_i (l + g^)(l + g^"+^) 

^ (1 - q^Y 


V 


0 


2 g- 


(1 -b g2"+4) 

(i-g 2)2 ; 


These results are g-analogues of some of the results given in §8.3]. Note moreover that 1T_(0) is a 
multiple of the identity, so that the commutant of W- equals the commuting algebra of Tirao m- Since 
the commutant is trivial, the weight W- is irreducible, which can also be checked directly. 


5. Quantum group related properties of spherical functions 

In this section we start with the proofs of the statements of Section |4| which can be obtained using 
the interpretation of matrix-valued spherical functions on Uq{g). 


5.1. Matrix-valued spherical functions on the quantum group. In this subsection we study some 
of the properties of the matrix-valued spherical functions which follow from the quantum group theoretic 
interpretation. In particular, we derive Theorem 14.31 from Remark 14.21 The precise identification with 
the literature and the standard Clebsch-Gordan coefficients is made in Appendix El and we use the 
intertwiner and the Clebsch-Gordan coefficients as presented there. 

We also need the matrix elements of the type 1 irreducible finite dimensional representations. Define 

^<3M2))-t C, = { Y { X ) ei , ei ,), n,m G 


where we take the inner product in the representation space for which the basis is orthonor¬ 

mal. Denoting 


t 


1/2 

- 1 / 2 ,- 1/2 
1/2 

'' 1 / 2 ,- 1/2 



7 


then a, /3, 7 , (5 generate a Hopf algebra, where the Hopf algebra structure is determined by duality of Hopf 
algebras. Moreover, it is a Hopf *-algebra with ^-structure defined by a* = 5, /]* = —gy, which we denote 
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by Aq{SU{2)). Then the Hopf ^-algebra Aq{SU{2)) is in duality as Hopf ^-algebras with Uq{su{2)). In 
particular, the matrix elements „ S Aq{SU(2)) can be expressed in terms of the generators and 
span Aq{SU{2)). Moreover, the matrix elements G Aq{SU{2)) form a basis for the underlying 
vector space of Aq{SU{2)). The left action of Uq{Q) on Aq{SU{2)) is given by {X ■ ^)(F) = ^(YX) 
for X,Y G Uqis) and ^ G Aq{SU{2)). Similarly the right action is given by • X){Y) = ^{XY) for 
X,Y G Uq{g) and ^ G Aq{SU{2)). A calculation gives „ and „ 

so that a ■ /3 • = q~^^'^l3, 7 • <5 • = q~^^'^S. 

Since and its powers are group-like elements of Uq{Q), it follows that the left and right action of 
and its powers are algebra homomorphisms. See e.g. m, m, m, m, and references given there. 
In the same way we view 

® tm2,n2 ' 

by taking the functions with respect to the tensor product Uq{Q) = Uq{su{2)) ®Uq{su.{2)). In particular, 
®tm 2 ,n 2 i^i^ 2 ) = ^mi .m <Jm 2 ,n 2 for A,/X G Similarly, the Hopf *-algebra spanned 

by all the matrix elements 0 is isomorphic to Aq{SU{2)) 0 Aq{SU{2)). We set Aq{G) = 

Aq{SU{2)) ® Aq{SU{2)) forG = SU{2)x SU{2). 

Define A = and let A be the commutative subalgebra of Uq{Q) generated by A and A~^. 

Recall the spherical function from Definition 14.41 and recall the transformation property (14.51) . 

Definition 5.1. The linear map $; Uq{Q) —^ End(H^) is a spherical function of type I if 

^{XZY) =t\X)^{Z)t‘^{Y), yX,YGB, yZGUgis). 

So the spherical function is a spherical function of type £ by (14.51) . 


Proposition 5.2. Fix £ G and assume (-^ 1 ,^ 2 ) satisfies (14.31) . 

(i) Write ® ^m,n! where the elementary matrices, then 




E 


E 


'~^mi,7712,722,n ^772i,ni 


2721,721 = —£1 7722,722 = —^2 


0 t 


A 

7722,722 ' 


(ii) A spherical function $ of type £ restricted to A is diagonal with respect to the basis {ep}p^_£. 
Moreover, for each A G Z, 

l.A(^^))m.n =E E 

i=-ti j=-l2 

so that £ 2 ( 1 ) is the identity. 

(iii) Assume $ is a spherical function of type £ and that 




t 

$= E $m.n0A^.„:G,(fl)^End(HO, 

with all linear maps onUq{Q) in the linear span of the matrix-elements t^il ®^i 212 ’ — *i’7i ^ 

£ 1 , —£2 < * 2 ,j 2 < £ 2 , then $ is o multiple 


Proof. Note (^>^,72 = (^ 1/2 for G Wg(fl), therefore we first compute 

Eor X G hlg(g) we have 


n 2 =— 


£1 


= (KjJ* E E ® 0 0 e: 

^7721 ,72i = —^1 7722,722= —^2 
£ £1 £2 

^^ 1 ,^ 2 ,^ ^^ 1 ,^ 2 ,^ ff£i 6d 

/ j / j 7721 ,7722 ,772'~^72i ,722 , 72 V*"7721,72i ^ ^7722,722 772 ’ 

772=—7721 ,72i = — ^l 7722 ,722 = —-^2 


/2 

7722 


This proves (i). 
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To obtain (ii) write $ = Ylln n=-t and observe = <^{A^Kf^) = 

^{A^)t^{Kf^) for all A e Z, /r e \TL implies since is diagonal. 

This gives ^m,n{A^) = 0 for m ^ n. Next pair with A^ using (i) and the observation made before 

Proposition [O] and the fact m = 0 unless nii — m 2 = m. Next take A = 0, and use (|A.4I) . 

Finally, for (iii) note that for X G ,B we have $(X) = t^(X)$(l) = $(l)t^(X). Since : B ^ End('H^) 
is an irreducible unitary representation, Schur’s Lemma implies that $(1) = cl is a multiple of the 
identity. Theorem 14.31 implies that for X G i3 


tHx) 


m,n 


E 


E 


^mi ,m2 ,m^ni ,n2 ,n V^mi ,ni 


1712,712 ——£2 




and, by the multiplicity free statement in Theorem l4.31 this is the only (up to a constant) possible linear 
combination of the matrix elements ® im 2 fixed {ii,i 2 ) which has this property. Hence, (iii) 

follows. □ 


The special case ^\/ 2 i/ 2 '- be calculated explicitly using the Clebsch-Gordan 

coefficients from (IA.5I) . Explicitly, 

V5= ^(9"^ +9)^?/2.i/2 

(5.1) = igP 1 _ 1 (g) 1 _ 1 - 1 ® 1 1 

Z 2’2 2’2 Z 2’2 2’2 

1 « 1 1 - 1 . r 

=- 9 a (g) a --/? (g) p --7 (g) 7 +-g 9 (g) d. 

This element is not self-adjoint in ^g(S'17(2)) (Si Aq{SU{2)). Recall the right action of Wg(g) on the map 
$: Uq{Q) —>■ End('H^), including the case 1! = 0, by ($ • X){Y) = $(Xy) for all X,Y G l^q(g)- This is 
analogous to the construction discussed at the beginning of this subsection. Then 

(5.2) -ip = ip- A~^ = ia (g) a - ® P — ^97 ^ J + ^ S = ip* 

is self-adjoint for the *-structure of 7l,j(5'17(2)) (S Aq{SU{2)). Then by construction 

(5.3) ip{iAXA-^)YZ) = e{X)iP{Y)s{Z), WX,Z G B, Y GUq{g). 

5.2. The recurrence relation for spherical functions of type £. The proof of Theorem 14.61 in 
the group case can be found in ESI Prop. 3.1], where the constants Aij are explicitly given in terms 
of Clebsch-Gordan coefficients. This proof can also be applied in this case, giving the coefficients Aij 
explicitly in terms of Clebsch-Gordan coefficients. A more general set-up can be found in [40l Proposition 
3.3.17]. The approach given here is related ESI Prop. 3.1], except that it differs in its way of establishing 
Xi/ 2 , 1/2 7^ 0. 

Proof of Theorem EJ As Uq ( 0 )-representations the tensor product decomposition 

j.i=±l/2 

follows from the standard tensor product decomposition for I4q{su{2)), see (lA.ll) . It follows that 

e 

E E 

±1/2 7n,n——£ 

where is in the span of the matrix elements G) ■ Note that 621), and a similar calcu¬ 
lation for multiplication by an element from B from the other side, show that has the required 

transformation behaviour (14. 5p for the action of B from the left and the right. Since the matrix elements 
4] sx®^r\ S 2 fo™ a basis for Aq{G), it follows that each 4'*’-^ satisfies (14.5E so that by ProT)osition l5.2l iiil 
Here = 0 in case {£i +i ,£2 + j) does not satisfy the conditions (14.31) . 

It remains to show that Ai/ 2 . 1/2 7 ^ 0. In order to do so we evaluate the identity of Theorem 14.61 at a 
suitable element of Uq{Q). For the Wq(su(2))-representations in (13.311 it is immediate that t^(e^) = 0 for 
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“ _ 1 G) _i + ^9 1 G) tf 1 

Z2’2 2’2 Z 2’2 2’2 







k > 2£ and that = 0 except for the case (r,s) = {—£,£) and then ^ 0. Extending to 

Uq{g) we find that = 0 if fci > 2ii or k 2 > 2 ^ 2 , and 




where the right hand side is non-zero in case m = —ii + £ 2 , n = £\ — £ 2 - 

So if we evaluate the identity of Theorem 14.61 at , it follows that only the term with 

ihj) = (1/2,1/2) on the right hand side of Theorem 14.61 is non-zero, and the specific matrix element is 
given by (El) with (£ 1 , £ 2 ) replaced by (£1 -I- £2 + 5 )- It sufhces to check that the left hand of Theorem 

14.61 is non-zero when evaluated at E 2 ^^^^. By (14.61) we need to calculate the comultiplication on 

E2ii+i e^^ 2 +i ^ Using the non-commutative g-binomial theorem m Exercise 1.35] twice we get 


A(^2^i-H1^2U-H1) ^ A{Ei)^E+1 + 1 


2ii + l2i2 + l 

E E 


■2£i-hl' 


■2£2-f r 

I 

_i 

g2 

-1 

CM 


^ k 2 — ^^ 2^2 + 1 —^2 


^ ^ 2 ^ 1 +l — fci ^ 2.^2 + ! —^2 


From (15.11) we find that (p{E^^E ^'^= 0 unless ki = k 2 = 1, and in that case the 
term /3 ® /3 of dSH) gives a non-zero contribution, and the other terms give zero. So we find 




■ 2 £i-hl' 


■ 2 £ 2 -£ 1 ' 

1 

,2 

1 


^{E,E 2 KIE ) (Ef ^^), 


and this is non-zero for the same matrix element (m, n) = {—£1 E £ 2 , £1 — £ 2 ) by (15.4E 


□ 


Note that we can derive the explicit value of Ai/ 2 . 1/2 from the proof of Theorem 14.61 bv keeping track 
of the constants involved. However, we don’t need the explicit value except in the case £ = 0. 

In the special case £ = 0, the recurrence of Theorem 14.61 has two terms and we obtain 


(5.5) 


— ^ 1 / 2 , 1/2 4>°j + i/ 2 ,^ 1 + 1/2 + ^- 1 / 2 ,- 1/2 4 >°j_i/ 2 ,^i_ 1 / 2 j 

1 1 - 
2®l-g4U-h2- 


, 1 _.l-q^E+4. 

^ 1 / 2 , 1/2 — X*? 


2^ l-qU^+2^ ^-l/2.-l/2 - 


The value of H. 1 / 2 , 1/2 can be obtained by evaluating at the group like element H'*', using Proposition l5.2f iil 
and 2tf{A^) = -|- q~^~^ and comparing leading coefficients of the Laurent polynomials in q^. This 

gives = ^i/ 2 ,i/ 2 (cEE/ 2 !E-i/ 2 .o)^’ ^ 1 / 2 , 1/2 follows from (IA. 6 E 

Having Hi/ 2 , 1 / 2 , we evaluate at 1, i.e. the case A = 0, which gives Ai/ 2 , 1/2 + ^- 1 / 2 ,- 1/2 = 5(9 + 9 ~^) 
by Proposition 15.2f iib from which A_if 2 .-i /2 follows. 


Proposition 5.3. For n G N we have 4)° 


l-g 2 


n,^n 4 ^ _ q2n+2 


UniF)- 


Recall that the Chebyshev polynomials of the second kind are orthogonal polynomials; 


(5.6) 


Un{cOs0) = 


sin((n -I- 1)0) 
shi 0 ’ 


Un{x)Ujn{x)\/^ - dx = dm.nE, 


xUji(^x'j — ^Uji^ii^x') -(- ^Un—l{x'j^ 


U-i{x)=0, Uo{x) = l, 


see e.g. [H], [ID], [DP] . 

Proof. From (15.51) it follows that $*] 1 = Pnip) for a polynomial Pn of degree n satisfying 

1 1 _ q2n+i 2 1 — g 2 ra 

XPn{x) = ^g~ \_^2 „+2 P"+l(^) + 

with initial conditions ^ 0 ( 2 ;) = 1 , Pi(a;) = q+q-i Set Pn{x) = g" i_qfL +‘2 rn{x), then ro(x) = 1 , 
ri(ai) = 2x and 2xr„(x) = rn+i(x) +rn-i(x). So r„(ai) = Un{x). □ 
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5.3. Orthogonality relations. In this subsection we prove Theorem 14.81 from the quantum group 
theoretic interpretation up to the calculation of certain explicit coefficients in the expansion of the 
entries of the weight. 

Recall the Haar functional on Aq{SU{2)). It is the unique left and right invariant positive functional 
ho : Aq{SU{2)) — 7 > C normalised by /io(l) = 1, see e.g. [9], [2TJ §4.3.2], [22, jlj. The Schur orthogonality 
relations state 


(5.7) 


ho: Aq{SU{2)) C, 




(1 - g«i+2) ■ 


Then the functional h = ho (S' ho is the Haar functional on Aq{G) = Aq{SU{2)) 0 Aq{SU{2)). We can 
identify the analogue of the algebra of bi-RT-invariant polynomials on G as the algebra generated by the 
self-adjoint element ijj, and give the analogue of the restriction of the invariant integration in Lemma 

EH 


Lemma 5.4. A functional r; Uq{Q) —>■ C satisfying T{fAXA~^)YZ) = e{X)'il){Y)e{Z) for all X,Z G B 
and all Y G Uq{Q) is a polynomial in if as in (15.21) . Moreover, the Haar functional on the ^-algebra 
C [■)/)] C Aq{G) is given by 


HpW) 



dx. 


Proof. From Proposition 15. 2l iiil and (15.3p we see that any functional on Uq{Q) satisfying the invariance 
property is a polynomial in ip, since this space is spanned by i ■A~^ for n G N. From Proposition l5.3l 

2 n, 2 71 

we have • A~^ is a multiple of Un{ip), since the right action of A~^ is an algebra homomorphism. 

The Schur orthogonality relations give = 0 for m n, and since 

^ Tl , ^ 71 ^ 771 ^ ^ 771 

the argument of h is polynomial in if, we see that it has to correspond to the orthogonality relations 
(EH) for the Chebyshev polynomials. So we find the expression for the Haar functional on the ^-algebra 
generated by '0- D 

Theorem 5.5. Assume h{q{g) —>■ End(lt^) are spherical functions of type £, see Definition 

Then the map 

t: Uqig) ^ C, X ^ tr((^ • • A-^)*){X), 

satisfies t{{AXA~^)YZ) = e[X)il){Y)e{Z) for all X,Z G B and all Y G h(q{g). 

In particular, any such trace is a polynomial in the generator if by Lemma l5.4l 
Corollary 5.6. Fix I G then for k,p G {0,1, • • • , 2i}, 

pAk 

W{ip)k,p ■■= tr(($^(o ^ ar{k,p) Uk+p-2rW, 

r—0 

with (lF(0)fe,p)* = W{ip)p^k- 

Proof of TheOrem\EE Write 4' = Y!m,n=-t '^rn,n ® = Ym,n=-i ® and 

linear functionals onUq{g), so that 


(5.8) 


e 

tr((T • H-i)($ • H-i)*) = ^ (^„,„ • 7l-i)($™,„ • H-i)* 

m,n——£ 

i 

T(y)= ^ ^vi.„,„(H-iy(i))4>™.„(H-i5(r(2))*), 

m,n——i (1^) 


using the standard notation A(y) = X)(y)^(i) ® ^( 2 ) and ^*{X) = ^{S{X)*) for Uq{g) —>■ C and 

XGUqig). 

Let Z G B and Y G h(q{g), then we find 


i 

t{YZ)= Y E '^mAA~^Y{i)Z(i))‘^mAA-^S{Y^2)rS{Z^2)r)- 

m,n=—t. {Y),(Z) 


Since H is a right coideal, we can assume that G B, so, by the transformation property (14.51) . 
'^m,niA-Ai)Z{i)) = T,i=-i'^rn,kA~^Y(^i))4AZii)). Using this, and 4_„(^(i)) = <,fe(^fi)) by the 
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^-invariance of B and the unitarity of and next move this to the <f>-part, and summing over n and the 
transformation property (14.51) for $, we obtain 

e i _ 

t{YZ) = E E E <^rnA^-^s{Y^2)rs{z^2))AiAAi)) 

m,k——l (y'),(Z) n——£. 

i _ 

m,k=-e. (Y) (Z) 

e 

= e{Z) Y. Y ci>fe,„(A-i5(y(2))*) = eiZMY), 

n,k=-t (Y) 

using X](z) '^(^(2))*^*!) = (S(z) ^{i)^A( 2 )))* = ^{Z) by the antipode axiom in a Hopf algebra. 

For the invariance property from the left, we proceed similarly using that A is a group-like element. 
So for Y G Uq{Q) and X G B we have 
i 

t(AXA-^Y)= Y E 'I'm.n(X(i)A-ly(i))$^,„(A-15(A)*^(X(2))*^(A-l)*,S(y(2))*). 

m,n——£ (X),(y') 

Now S'(A)* = A~^, S'(A“^)* = A. Proceeding as in the previous paragraph using X(i) G B we obtain 

i s. _ 

t{AXA-^Y) = E E E 

n,k=-i {X),{Y) m=-i 

I 

= E E Y ^kAY^)A-^s(XmyAS{Y^-,^Y). 

n,k=-t. (Y) (X) 

The result follows if we prove X](x) ^(*i)^~^‘^(X2))*^ = A~^£(X). In order to prove this we need the 
observation that S(X*) = A~‘^S{X)*A^ for all X G Uq{Q), which can be verified on the generators and 
follows since the operators are antilinear homomorphisms, see Remark (EZl). Now we obtain the required 
identity; 

= = e{X*)A-^ = 7{X)A-\ □ 

(^) (X) 


Remark 5.7. The required identity S{X*) = A~'^S{X)* for all X G Uq{Q) can be generalised to 
arbitrary semisimple g. Indeed, since the square of the antipode S is given by conjugation with an 
explicit element of the Cartan subalgebra associated to p = ^ Sa>o ®-S- ESI Ex. 4.1.1], and since 

in a Hopf ^-algebra S' o * is an involution, we find that S{X*) = K-pS{X)*Kp if S‘^{X) = K-pXKp as 
in |331 Ex. 4.1.1]. 


n,*n ' 

0 . A-1 

n, in 


Proof of Corollary 15.61 By Theorem l5.5l and Lemma [5.4l the trace is a linear combination of 

2 

hence a polynomial in ip. To obtain the expression we need to find those n’s for which i ■ A 
occurs in W{ip)k,p by Proposition 15.31 Using (14.4L (15.811 and Proposition 15.21 we find that in W{ip)k,p 

—k —V £——k £——v 

matrix only matrix elements of the form di )* ® ^a 2 fc 2 ^2 Using the Clebsch-Gordan 

decomposition we see that bii^ci di)* similarly )* can be written as a sum of matrix 

elements from j. g {0,1,--- , fcAp}, and similarly g g {0,1,-- - ,{2£—k)A(2i—p)}. 

By Proposition l5.31 and Proposition l5.21 the only n’s that can occur are n = k+p—2r, r G {0,1, • • • , fcAp}. 
The statement on the adjoint follows immediately from (15.81) and ip being self-adjoint, see (15.21) . □ 


We now can start the first part of the proof of Theorem 14.81 except for the fact that we have to 
determine certain constants. This is contained in Lemma [ESI 


Lemma 5.8. IFe have 


E E E ^ q 


(i-g2) ■ 


i——£ m\ ——£\ 1712 ——£ 

The proof of Lemma 15.81 is a calculation using Theorem 14.61 which we postpone to Subsection I6.3[ 
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First part of the proof of Theorem \4-ttl Using the notation of Section 0] we find from Corollary 14.71 and 

(BH) 


21 2i 

k—Q p—0 

2& 21 

= EE 

k=0p=0 


A-y) 



where we use that the action by >1 ^ from the right is an algebra homomorphism, since A ^ is group 
like, and that if is self-adjoint. By Proposition [53 Lemma [5.41 and (EH) we have 


(5.9) 


J-1 

= 5, 


[{P„{x)rWix)Pmix))^yi-x^dx = h(tr((ci>f(„,,) . 


g2^+2"(l _ g2^2 


£ 


^ Q2n-\-2i-\-2'^ _ q4:£-\-2n — 24 - 1 - 2 ^ 

£+^(n—i) 


i(n-l-i) 

E E 


E .„^{n+i),£+^(n-i),l,2 2(a+b) 

\'-'a,b,r I 9 



after a straightforward calculation. Plugging in Lemma l5.8l in (15.9L and rewriting proves the result using 
Corollary [53 □ 


Note that we have not yet determined the explicit values of a*( to, n) in Theorem l4.8l and we have not 
shown that lU is a matrix-valued weight function in the sense of Section |2j The values of the constants 
at{m,n) will be determined in Section o and the positivity of lU(a;) for x S (— 1 , 1 ) will follow from 
Theorem 14.151 


5.4. (/-Difference equations. It is well-known, see e.g. Koornwinder m, Letzter [3S], Noumi EZ], 
that the centre of the quantised enveloping algebra can be used to determine a commuting family of q- 
difference operators to which the corresponding spherical functions are eigenfunctions. In this subsection 
we derive the matrix-valued g-difference operators corresponding to central elements to which we find 
matrix-valued eigenfunctions. 

The centre of Wq(g) is generated by two Casimir elements, see Sectional 


Ui = 


qKf 


^-1 


Kt -2 


(g-g-i)2 

Because of Proposition E3 and (13.41) we find 


-I- Ui Fi, 


U 2 = 


qKf 


^K2-2 


{q-q-^y 


-|- E 2 F 2 . 


(5.10) 


U,; • $1 » = $ 


£ 1,(2 


n, = 



i = 1 , 2 . 


The goal is to compute the radial parts of the Casimir elements acting on arbitrary spherical functions 
of type £ in terms of an explicit g-difference operator. In order to derive such a g-difference operator, 
we find a SylS-decomposition for suitable elements in Uq{g) in Proposition 15.101 This special case of 
a ,B.4,B-decomposition is the analogue of the ifAiL-decomposition, which has not a general quantum 
algebra analogue. For this purpose, we first establish Lemma 15.91 which can be viewed as a quantum 
analogue of [ 8 j Lemma 2.2], and gives the SylS-decomposition of 

Lemma 5.9. Recall A = . For A G Z \ {0} we have 

P^A^ = {k^/'^A^Bi - q^RP^BiA^^ . 

Proof. Recall Definition 14.11 so the result follows from 

A^Bi = + qA^E2Kf^''^K]!'^ = g(g^“^ - (y"^^)Kf^l'^Ky'^E2A^ F q^B^A^, 

and using e B. □ 
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Proposition 5.10. The Casimir elements Hi and H 2 admit the BAB -decompositions 

^ ^ g(l - Xl/2 /^A+l _ 4 A , g3(l-gr2A) ^A-1 

^ (l-q2)2(l_g2A+2) (l-g2)2 ^ (1 _ ^2)2 (1 _ ^2A+2) 


« 2 A +1 


(1 - g2A+2)2 

+ (l_g 2 A+ 2 ) 2^1 


« 2 A +2 


(l_q2A+2)2- 


„3A+3 

(1 — q 2 A+ 2^2 


and 


noA^ = 


g(l_^2A+4) 2g2 


A 'd\^ 'J ) „_i/2 .A+1 _ ^g aA , g^(^ g^^) r^l/2 ^A-l 

(l-q2)2(l_g2A+2) (l-<?2)2 (1 _ g2)2 _ ^2A+2) 


„2A+1 «2A+2 

® jB2K^^^BiA^+^ - —^ 


(1 - g2A+2)2 


A^+^K^ABiB2 


(l-g 2 A+ 2 ) 2 ^ 

^A „3A+3 

Proof. We first concentrate on H 2 , the statement for Hi follows by flipping the order using cr as in 
R,emark l4.2f iv'l. We need to rewrite E 2 F 2 A^, which we do in terms of F 2 A^ and next using Lemma IhTiH 
The details are as follows. 

Using Definition 14.11 the commutation relations and pulling through F 2 to the left, we get 
(5.11) B2F2A^ = q-^F2F2A^-^ + q^FiF2K^^‘^A^. 

Similarly, and only slightly more involved, we obtain 


(5.12) 


F2A^B2 = q^-^F2F2A^-^ - q ^-^——+ q^-^FiF2K^/‘^A^. 


q-q 


Using (15.111) and (15.121) we eliminate the term with F 1 F 2 , and shifting A to A + 1 gives 

(5.13) {q-^ - q^^+^)F2F2A^ = ^2^224^+^ - q^+^F2A^+^B2 - g2A+l -^2 - ^A_ 

q-q 

Apply Lemma 15.91 on the first two terms on the right hand side of (15.131) and note that the remaining 
terms in (15.131) and in HzA'*' can be dealt with by observing that K 2 = K~^I‘^A = AK~^I'^. Taking 
corresponding terms together proves the SA-B-decomposition for HzA^ after a short calculation. □ 


Our next task is to translate Proposition (STTU] into an operator for spherical functions of type £, from 
which we derive eventually, see Theorem l4.131 an Askey-Wilson g-difference type operator for the matrix¬ 
valued orthogonal polynomials P„. Let $ be a spherical function of type 1., then we immediately obtain 
from Proposition 15.1()1 

(5.14) 

g3(l-g2A) 

^ (l-(Z2)2(l-g2A+2)M^ ) 

9A-L9 

.$(2l^+l)i^(ii:-l/2^2Pi) 


„2A-|-1 o2A-|-2 

^ t^(PiiWi/2^2)$(A^+i) - ^ 


(1 _ g2A-|-2)2‘' ^ 2 ; V ; ^ 2 ^ _ ^2A-|-2)2 

„A „3A+3 


(1 - g 2 A+ 2 ) 


The analogous expression for ^(HzA^) of (15.141) can be obtained using the flip a, see Remark I4.2f ivl. 
In particular, it suffices to replace all t^{X) in (I5.14p by J^t^{X)J^, see Remark 14.21 ^1. to get the 
corresponding expression. 

By Proposition I5.2r iil we know that $(A^) is diagonal. We can calculate the matrix entries of 
$(HiA'^) using the diagonal matrices t^{K^^^'^), t^{BiK~^^'^B 2 ), t^{K~^/'^B 2 B 1 ) and the the upper 
triangular matrices t^{BiK~^/'^), t^{Bi) having only non-zero entries on the superdiagonal and the lower 
triangular matrices t^{K~^^'^B 2 ), ^^(^ 2 ) having only non-zero entries on the subdiagonal, see (021), dSH). 
By (I5.14|) . $(HiA^) is a diagonal matrix as well. 
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For $ a spherical function of type I we view the diagonal restricted to ^ as a vector-valued function 






m——£ 

So we can regard the Casimir elements as acting on the vector-valued function and the action of the 
Casimir is made explicit in Proposition 15.Ill 

Proposition 5.11. Let ^ be a spherical function of type I, with corresponding vector-valued function 
$: —>■ representing the diagonal when restricted to A. Then 

(1 - q^y 

where Mf^z) is a tridiagonal and ^^ 1 ( 2 :) is a diagonal matrix with respect to the basis of TL^ 

with coefficients 

(aA( M q^+^il-zy 


{Mi[z))m,m+1 — 
iM({z))m,m = 


(l-g 2 ) 2 (l_g 2 ^ 2 )’ 

-„m+l 

ybym + i)y, 


(1 - q^z^y 


^2q2+m 


(l-g2)2(l_g2^2) (1_ ^2^2)2 

2 


{{hym)y + {bym + i)y) 


^3^m+3 

{Myz))m,m -1 = \ ,,J b\m)y. 


y-q^z^y' 


Moreover, for LI 2 the action is 

$( 112 ^^) = Mi(qy^{A^+y - + Ni{qy^{A^-y, 

(1 - g 2)2 

where M^z) = J^Mi{z)J^ is a tridiagonal, and M^z) = J^M({z)J^ is a diagonal matrix. 

Remark 5.12. The proof of Theorem 14.1.11 does not explain why the matrix coefficients of pq and Pq-i 
in Theorem 14. Ill are related by z -O- z~^. In Proposition [STT] there is a lack of symmetry between the up 
and down shift in A, and only after suitable multiplication with $0 from the left and right the symmetry 
of Theorem 14.111 pops up. It would be desirable to have an explanation of this symmetry from the 
quantum group theoretic interpretation. Note that the symmetry can be translated to the requirement 
4'(z)iVi(g-2z-i) = Mi(z)4'(gz) with 

The remark following (I5.14p on how to switch to the second Casimir operator gives the conjugation 
between Mf and M|, respectively Nf and 1V|. Note that in case £ = 0, $ and $ are equal, and we find 
that Proposition [531] gives the operator 

= ^{niAy 


-- 

(l-g2)2(l_^2A-h2) V ) (l-q2)2 V 


g3(l-g2A) 


.$(21^-4) 


(l-g2)2(l_^2A-h2) 

for $: Uq{Q) —>■ C a spherical function (of type 0), which should be compared to (SO] Lemma 5.1], see 
also [3Sj, [57] . 

Proof. Consider (15.141) and calculate the (TO,m)-entry. Using the explicit expressions for the elements 
t^{X) for X G B, see (14.21) . (II.ip . in (I5.14p . we find 

., 2 A^ 0„2 


.p.o .An _ .J./.A- 1 N V 

)m,m — JZ - 0 . 2/1 - 2A-I-2N )m,m — 


(l-g2)2(l_g2A-H2) 


(1 - 


+ 




(l-g2)2(l-g2A-h2) 
2A-|-2-|-77i / ^ ^ 2 


(1 - g 2 A-H 2)2 

f l,£. f rr-y-t I 1^^2 


<?2A+2+™(6^(m))2 q^+-^+yb\m + l)y 

(1 _ g2A-|-2)2 )m,m+ - g2A-|-2)2 


q' 


3A-|-3-|-?ti / 1-.£ 


{bym)f 


(1 


_ „2A-|-2'|2 
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which we can rewrite as stated with the matrices M({q^) and N({q^). The case for ^2 follows from the 
observed symmetry, or it can be obtained by an analogous computation. □ 


In particular, we can apply Proposition 15.111 to using (jb.lOll to find an eigenvalue equation. 

In order to find an eigenvalue equation for the matrix-valued polynomials, we first introduce the full 
spherical functions 

u 

(5.15) $1: ^ ^ End(C2^+i), ^ ® E.,,, 

z,j=0 

for n S N. So we put the vectors o)> ‘‘ it) columns in a matrix, and we relabel in order 

to have the matrix entries labeled by i, j S {0, • • • , 2£}. We reformulate Proposition 15.Ill and (15.1011 as 
the eigenvalue equations 




(5.16) 


2i 


A.( 1 ) = E 


l-n-j I „ 3 +n+j 


i=o 


(1 - 


-Ejj, A„(2) = JA„(1)J, 


where (Mi{z))^ ^ for m,n £ {0,1, • • • , 2i}, and similarly for Ni, are the matrices of 

Proposition 15.111 shifted to the standard matrix with respect to the standard basis {e„}^^Q of 
Note that the symmetry of Proposition [5TT] then rewrites as M 2 {z) = JMi{z)J, N 2 {z) = JNi{z)J, with 
J' tn I y 62 £—n- 

Now we rewrite Corollary 14 . 71 after pairing with as 


(5.17) 


= ¥^{AYPn{^i{q^+)) £ End(C 2 '+i), 


where /i(a:) = ^{x + x ^), using that pairing with the group-like elements A^ is a homomorphism and 
ip{A^) = /r(g^+^) by (15.111 . Using (I5.17P in ()5.16p proves Corollarv l5.13l 

Corollary 5.13. Assuming ^q(A^) is invertible, the matrix-valued polynomials Pn satisfy the eigenvalue 
equations 

Pnitiiq^)) A„(z) = MYq^) PniKq^^Y) + N^{q^) Pn{li{q^-)), ^ = 1, 2, 


MYq^) = {¥,iA^-Y) ^M,(g>-i)cI>^o(A>), NYq^) = ($g(A^-i)) ^iV,(q^-i)<i>g(A>-2), 
where A„(f) are defined in (15.1611 and p,{x) = ^{x -\- x~Y- 


It remains to prove the assumption in Corollary 15.131 for sufficiently many A, and to calculate the 
coefficients in the eigenvalue equations explicitly. This is done in Section [71 

Having established (I5.17L we can prove Corollary I4.9l bv considering coefficients in the Laurent ex¬ 
pansion. 

Proof of Corollary \4.9[ The left hand side of (15.1711 can be expanded as a Laurent series in q^ by Propo¬ 
sition I5.2f iil and (I5.15P . The leading coefficient of degree n is an antidiagonal matrix 

(41(^5).,, = + lower order term., 

since the Clebsch-Gordan coefficient is zero unless i-£ j = 21. With a similar expression for $q(A''') on 
the right hand side and expanding Pn{lJ'iq^~^^)) = lc(P„)g”^^+^)2“" -|- lower order terms gives 

-,h{n+j),e+i(n-j),e 


]„dp 'I _ r no —rt /pt 2 2 (” \ f \ 


and (IA.7II gives the result. 

Corollary 14.91 gives Po{x) = /, so Corollary 15.131 gives 

Ao{l)=MfiqY+NfiqY, VAeZ. 


□ 


(5.18) 
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5.5. Symmetries. Even though we can use the explicit expression of the weight W to establish Propo¬ 
sition UTTUl we show the occurrence of J in the commutant from Remark 14. 5f iil. 

Observe that (£i,£ 2 ) = ^(n,k) gives (^ 2 ,^ 1 ) = ^{n,2i — k) and that (t{A) = A, so Remark [d.bi iil 
yields • 4l“^)(cr(^)) = {^l(n, 2 e-i) ‘ A-^){Z)J^ for any Z S Ugis). Similarly, using moreover 

that (T is a ^-isomorphism and that S and a commute, see (1331) we obtain • A ^)*{a{Z)) = 

2 i-i) ■ A~^Y {Z)J^. Since {a ® a) o IS. — a and (J^)^ = 1, we find for Z G Uq{Q) from these 
observations, cf. (1531) . 

By the first part of the proof of Theorem 14.81 

ja(Z)) = ((P„(V^))*iE(V^)P™(V^))^^_^^^_ jZ). 

Note that ip{a{Z)) = ip{Z) by the symmetric expression of (15.21) . Again using (cr 0 cr) o A = A o cr, we 
have p{ilj){cr{Z)) = p{'>p){Z) for any polynomial p. It follows that 

(5.19) ((P„(V^))*IE(V^)P™(V^))^^ = 

For n = m = 0, (15.191) proves that J is in the commutant algebra for W as stated in Proposition 14. 101 
Note that (15.191) . after applying the Haar functional on the ^-algebra generated by ^ as in Lemma 
1331 also gives JG„J = G„ as is immediately clear from the explicit expression for the squared norm 
matrix G„ in Theorem 14.81 derived in Section |5.31 It moreover shows that (JP„J)„gN is a family of 
matrix-valued orthogonal polynomials with respect to the weight W. It is a consequence of Corollary 
14.91 see Section[31 that JP„J = Pn, since Jlc(P„)J = lc(P„) by Corollary 14.91 

Note that we have now proved Proposition 14.101 except for the D-inclusion in the first line. This is 
done after the proof of Theorem 14.81 is completed at the end of Section 16.11 

As a consequence of the discussion on symmetries, we can formulate the symmetry for in Lemma 

ra 

Lemma 5.14. With J: en^ e 2 i-n we have J = for all A G Z. 

Proof. This is a consequence of initial observations in Section 331 For i,j G {O,--- ,2£}, using (15.151) 
and <j{AY = A^ we obtain 

(J4'„(yl")j). , = = {•S-i,,.,,,-!,(A’-)), 

6. The weight and orthogonality relations for the matrix-valued polynomials 


In this section we complement the quantum group theoretic proofs of Section [5] of some of the state¬ 
ments of Section 2] using mainly analytic techniques. In Section 16.11 we prove the statement on the 
expansion of the entries of the weight function in terms of Chebyshev polynomials of Theorem 14.81 In 
Section 33] we prove the LDU-decomposition of Theorem l4.15l In Section 3751 we prove Lemma 3751 using 
a special case and induction using Theorem 14.61 in the induction step. 


6.1. Explicit expressions of the weight. In order to prove the explicit expansion of the matrix entries 
of the weight W in terms of Chebyshev polynomials, we start with the expression of Corollary 15.61 for 
the matrix entries of the weight W. After pairing with A'*', we expand as a Laurent polynomial in in 
Proposition 16.II Then we can use Lemma 16.21 whose proof is presented in Appendix IB.II 


Proposition 6.1. For 0 < k,p < 2£, X G Z we have 

i(fe-l-p) 

W{'if)k,piAY = di{k,p) = di^{k,p), 

s=-^{k+p) 

ifc ip e 

^ ^ ^ ^2(z-|-j-n)+2(i-|-ife)(i-n+^-ifc)+20 + ip)0-n-|-£-ip) 

i=-^k j=-^pn=-t. 
s=j-i 


k 


to 

1 


P 


to 

1 


2£ 

\k-i 

9 " 

£ — ^k + n — i 


Ap-j. 

9 " 

_£- \p + n-j 

,2 

£ — n 
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Proof. We obtain from Corollarv l5.6l 


WUh.M^) = j;: 


(6^1) 


m,n——t 




using that is a group like element and 5'(A^)* = A Proposition l5.2f ii') gives the explicit expression 




= E E E g2A0-d^2(,+,-n)^ 

n=-e i=-lk j=-^p 

where the Clebsch-Gordan coefficients are to be taken as zero in case \i — n\ > £ — respectively 
\j — n\ > i — ip. Now put s = j — i, then s runs from —^{k+p) up to ^{k+p), and we have the Laurent 
expansion 


i(fc+p) 

WWkA^^) = E ^'sik,p)q 

S=-i(fc+p) 


2sX 


with 


2P e 


<iKp) = 1 : 1 : 1 : ( 4 “:’“) ( 44 ;!"') 

i=-^k j=-^pn=-t 


Plugging in the explicit expression (IA.3I) gives the explicit expression for dl{k,p). 

In order to show that df (p, k) = df_g{p, k), we note that p{'ip){A^) = p{p{q^)) is symmetric in A O —A 
for any polynomial p and so Corollary 15.61 implies the symmetry. □ 


Note that a similar expression can be given for a matrix element of P„(4)*IP('4)4m('4)) cf. the first 
part of the proof of Theorem 14.81 but this is not required. The symmetry in the Laurent expansion of 
W{'ij:)k,p{A^) does not seem to follow directly from known symmetries for the Clebsch-Gordan coefficients, 
see e.g. [2T1 Ch. 3]. 

Now we can proceed with the proof of Theorem 14.81 for which we need Lemma 16.21 

Lemma 6.2. For £ S and for k,p G {0, • • • , 2£} subject to k < p and k + p < 2£ we have with the 
notation of Proposition 16.11 and Theorem \4.8[ 

k k-\-p—2r 

E E o:r{k,p) = dl{k,p). 

r—0 a—0 
2s—k-\-p—2{r-\-a) 

Lemma 16.21 contains the essential equality for the proof of the explicit expression of the coefficients 
of the matrix entries of the weight of Theorem 14.81 The proof of Lemma 16.21 can be found in Appendix 
m and it is based on a g-analogue of the corresponding statement for the classical case [531 Thm. 5.4]. 
Previously in 2011, Mizan Rahman (private correspondence) has informed one of us that he has obtained 
a g-analogue of the underlying summation formula for m Thm. 5.4]. It is remarkable that Rahman’s 
g-analogue is different from the one needed here in Lemma 16.21 


Second part of the proof of Theorem |4.d[ We prove the statement on the explicit expression of the matrix 
entries of the weight in terms of Chebyshev polynomials. By Corollary 15.61 we have 

k/\p kAp k-\-p—2r 

W{il))k,p{A^) = ^ar{k,p)Uk+p-2T{p{q^)) = ^ E 

r=0 r—0 a—0 


i(fc-l-p) 

= E 

S = -i(fe+p) 




kAp k-\-p—2r 

E E <^Ak,p) 

r—O a—0 
\2s—k-\-p—2{r-\-a) / 


^2sA 
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for all A S Z. Since the coefhcients ar{k,p) are completely determined by W{'ip)k,p and since W{il^)k,p = 
2£-k = {w{^)p,ky , we can restrict to the case k < p, k + p < 2£. For this case the result 
follows from Lemma 16.21 and hence the explicit expression for ar{k,p) in Theorem l4.8l is obtained. □ 

Note that proof of Theorem 14.81 is not yet complete, since we have to show that the weight is strictly 
positive definite for almost all a:: S [—1,1], see Section [5J This will follow from the LDU-decomposition 
for the weight as observed in Corollary 14.161 but in order to prove the LDU-decomposition of Theorem 
14.151 we need the explicit expression for the coefficients at{m,n) of Theorem 14.81 

In Section [nm we observed that J commutes with W{x) for all x. In order to prove Proposition 14. 101 
we need to show that the commutant is not larger, and for this we need the explicit expression of at{m, n) 
of Theorem l4.8l 

Proof of Proposition \4-10\ Let Y be in the commutant, and write W{x) = ^kUn{x)^ with Wk G 

Mat 2 ^+i(C) using Theorem 14.81 Then [Y, lUfc] = 0 for all k. The proof follows closely the proofs of 
ESI Prop. 5.5] and ESI Prop. 2.6]. Note that W 2 t and W 2 t-i are symmetric and persymmetric (i.e. 
commute with J). Moreover {W 2 t.)m,n is non-zero only for m n = 2£ and (W 2 £-i)m,n is non-zero 
only for \m + n — 2£\ = 1. From the explicit expression of the coefficients at{m,n) we find that all 
non-zero entries of W 2 i and W 2 ^_i are different apart from the symmetry and persymmetry. The proof 
of Proposition 14. 101 can then be finished following [231 Prop. 5.5]. □ 


6.2. LDU-Decomposition. In order to prove the LDU-decomposition of Theorem 14.151 for the weight 
we need to prove the matrix identity termwise. So we are required to show that 

mAn 

( 6 . 2 ) ^ ^ 

for the expression of W{x) in Theorem 14.81 and for the expressions of L{x) and T(x) in Theorem 14.151 
Because of symmetry we can assume without loss of generality that m > n. Then (16.21) is equivalent 
to Proposition 16.31 after taking into account the coefficients in the LDU-decomposition, so it suffices to 
prove Proposition 16 . 31 in order to obtain Theorem 14.151 


Proposition 6.3. For 0 < n < m < 21, I S and with at{m,n) defined in Theorem \4-8[ we have 

n n / ^ 2fc+2 j 2\ 

atjm, n)Um+n- 2 tix) = y^/3fc(m,n) ^ 




(3kim,n) 


iq^;q^)n 2 ,,^,y-,q^)2i+k+liq^;q^)2i+k 

{q^-,q^)m+k+l{q^;q^)n+k+y'^ ^ q^^{q^;q^)l, 


Before embarking on the proof of Proposition 16.31 note that each summand on the right hand side of 
the expression of Proposition 16.31 is an even, respectively odd, polynomial for m + n even, respectively 
odd, since the continuous 9 -ultraspherical polynomials are symmetric and since 


k 


i=i 

see (|4.12|) . is an even polynomial with a factor (1 — x^). In the proof of Proposition 16.31 we use Lemma 

O 


Lemma 6.4. Let 0 < k < m < n and t gN, then the integral 
1 f w{x q 
J_i vl — 

is equal to zero for t > m and for 0 <t < m it equals Ck{ra, n)Rk{fJ'{t); 1,1, q~‘^™'~^, q^) with 

q-2{1) (^ 2 fc+ 2 . (g 2 fcH- 2 . (^ 2 ^+ 2 . ^ 2 )_^ 


Ck{m,n) = 


1 _ q2k+2 (q,4fc-i-4. {q^-,q^)m-k {q^-,q^)n-k 

(-l)^(g^*^+^; q^)k+i 

{q^\q^)m+n+l 
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In Lemma 16.41 we use the notation (14.131) for the g-Racah polynomials. Lemma 16.41 shows that the 
expansion as in Proposition 16 . 31 is indeed valid, and it remains to determine the coefficients /3fc(m,n). 

The proof of Lemma EH] follows the lines of the proof of [241 Lem. 2.7], see Appendix IB.21 The main 
ingredients are, cf. the proof in [24], the connection and linearisation coefficients for the continuous 
g-ultraspherical polynomials dating back to the work of L.J. Rogers (1894-5), see e.g. [H §10.11], [20l 
§13.3], [ini (7.6.14), (8.5.1)]. So we write the product as a sum over 

r of continuous g-ultraspherical polynomials Cr{x; using the linearisation formula and we write 

Un+m.- 2 t, which is a continuous g-ultraspherical polynomial for /3 = g, in terms of C's(a;;g^^“''^jg^) using 
the connection formula. The orthogonality relations for the continuous g-ultraspherical polynomials then 
give the integral in terms of a single series. The details are in Appendix lB.2l From this sketch of proof it 
is immediately clear that Lemma ElTl can be generalised to a more general statement. This is the content 
of Remark 1 6. 51 whose proof is left to the reader. 


Remark 6.5. For integers 0<t, k<m<n and parameters a, /3, we have 


1 r w{x q q'^+^a\q)Cjn+n- 2 t{x; ^\q)dx 

27r .^1 - 

= Cfc(TO,n, a,^)4<P3 


where 


Ck{m,n,a,l3) = 


a ^q ^,g ^,ag^+^,/3a ^g ‘ ^ 

a“ig“™, a“ig“”,/3 

(a2g2fc+2, g, gfc-‘+l, a-2g-™-"-l, /3-lg'=-™-"+‘, /3-lg; q)^ 

^ (Qg^+^;g)m-fc {aq'^+^;q)n-k (a^g^^+^;g)m-en-2fc 

{q;q)m-k {q;q)n-k {aq>^+‘^; q)m+n-2k 

ip- _ \Q)k-t iP]q)m+n—t-k , h+l \k-t 

{q;q)k-t {aq'^+'^; q)m+n-t-k^ 


Note that the 4 (^ 3 -series in Remark 16.51 is balanced, but in general is not a g-Racah polynomial. 

In the proof of Proposition 16.31 and hence of Theorem 14. 1 51 we need the summation formula involving 
g-Racah polynomials stated in Lemma 16.61 Its proof is also given in Appendix IB. 21 


Lemma 6.6. For I S and m,n, k G N with 0<k<n<m we have 

(.,2m-U.,,2\ (^4i+4-2t. ^2\ ,,, 






(g2™+4;g2)„_( {q^]q^)t 


= <7 


i{ri — l) — k(k+l)—4ni / _i \n-t-A: (Q )2£+k+l{Q :Q ^2£—k _ (^ Q _)_ 

^ ^ iq^;q^)2e+l {q^; q^Uq^; q^hi-n 


With these preparations we can prove Proposition l6.31 and hence the LDU-decomposition of Theorem 

m 


Proof of Proposition 16.51 Since we have the existence of the expression in Proposition 16.31 it suffices to 
calculate /3fe(m,n) having the explicit value of the at(m,n)’s from Theorem 14.81 Multiply both sides 
by — x‘^Um+n- 2 tix) and integrate using the orthogonality for the Chebyshev polynomials, so that 

Lemma 16.41 gives 

1 ” 

-atim, n) = ^ /3fc(m, n)C'fc(m, n)Rk{^iit); 1,1, g“^”"“^, g^). 

k^O 

The g-Racah polynomials Rfc(//(t); 1,1, satisfy the orthogonality relations 


^-2^-2^ ^-2n-2. 1, 1, q^) 




/-| „2m-|-2\/-| „2n-|-2 /„2m-l-4 „2n-|-4, „2\ 

^ r - 2 /c(m+n+i) (1 - g ^ )(i - g (g ^ ,q ^ ,q )k 

(l_g4fe+2) ^q-2m^q-2n.q2)^ ' 
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Using the orthogonality relations we find the explicit expression of /3fe(m,n) in terms of at{rn,n), and 
using the explicit expression of at{m,n) of Theorem 14.81 gives 


^ q 2 k{m+n+l) _ ^ 4 fe+ 2 ^ 

I3k{m,n) = - 


2m 2n.„2'\ 

[Q ,Q ,Q )k 


4 Ck{m,n) (1 — g2™+2)(l — g2n+2^ ('^2m+4^ ^2n+4. 

q 


^24^2^ _ ^2m+2n+2-4t 


)Rk{Kty, 1,1, q-^^-y g-"”-"; q^)at{m,n). 


t=o 


This expression is summable by Lemma 16.61 Collecting the coefficients proves the proposition. 


□ 


Last part of the proof of Theorem \4-8\ Now that we have proved Theorem l4.151 Corollary 14. 161 is imme¬ 
diate, since the coefficients of the diagonal matrix T{x) are positive on (—1,1). So the weight is strictly 
positive definite on (—1,1), which is the last step to be taken in the proof of Theorem 14.81 □ 

6.3. Summation formula for Clebsch-Gordan coefficients. In this subsection we prove Lemma 
15.81 which has been used in the first part of the proof of Theorem 14.81 see Section 15.31 The proof of 
Lemma 15.81 is somewhat involved, since we employ an indirect way using induction and Theorem 14.61 


Proof of Lemma \5.H\ Assume for the moment that 


(6.3) 




2 i _ 2 ii- 2 i- 2 m 


q = q 


(1 - g4^+2) 

_ n4<!i-|-2 


(1 


)■ 


Assuming (16.31) the lemma follows, using ^ = 0 if mi — m 2 f. 


E \ ^ \ ^ |2 2 (mi+m 2 ) _ \ ^ „ 4 mi \ ^ 




= E 


q 


2 ei- 2 i+ 2 m 




. 4 :£-\- 2 \ 


(1 - g4U-H2) 


= q 


(l-g2) ■ 


It remains to prove (16.31) . We do this in case -b ^2 = Put {£ 1 , 12 ) = {k/2,£— k/2) = ^{0,k) for 
k G N and k < 2£. Using the explicit expression for the Clebsch-Gordan coefficients (1A.31) . we find that 
in this special case the left hand side of (ESI) equals 

y-_ {q'^;q'^)k{q'^;q^)2e-k{q'^;q^)i-^iq^;q'^)e+i _ ^.2i+2(m-kl2)(m-i+e-kl2) 

^ isP'\q^)ki2-m{q^\q'^)ki2+miq'^;q‘^)e-ki2-m+iiq‘^]q‘^)i-ki2+m-i{q^] q^)2i 

where the sum runs over i for —£ + k/2 + m<i<£—k/2 + m. Substitute i>-^p — £ + k/2 + m to see 
that this equals 

i—k/ 2 -\-m 

E 


{q ]q )£-iiq ',q )e+i 


— 2 i{l-\-m-\-k/ 2 ) 


i——i-\-k/2-\-m 


{fP" 1 P')t—k/2^m—pPi Q‘^')i—kj2—m-\-i 


E jq ) 2 e-p-k/ 2 -m{q )p+k/ 2 +m 2 (l+m+k/ 2 )(p-e+k/ 2 +m) 

{qyqW:Q^)2e-p-k ^ 


Simplifying the expression we find 

— 2 (l+ 77 i+fc/ 2 )(—£+fc/ 2 +m) 2 £—k /2 — miQ -iQ )fc/ 2 +'r 


{q‘^-q^)2i-k 


^—Al-\- 2 k „ 2 m+fc +2 

, y 2 - 2 k -2 

■2Pl\ ^.u+k+2m ,g 


_ - 2 (l+m+k/ 2 )(-e+k/ 2 +m) ^^ ’g ) 2 £-fc/ 2 -?n(g (9 )k/ 2 +Tn jq ‘^]q^) 2 l-k 


= q 


{q^;q^)2i-k (g-4^+'=+2™;Q2)2^-fe’ 

since the 2^1 can be summed by the reversed g-Chu-Vandermonde sum |161 (11.7)]. Putting everything 
together proves (I01) . and hence Lemmafor the case {£i,£ 2 ) = C(0, k), i.e. £1 + £2 = £■ 

To prove Lemma 15.81 in general, we set 




c: 


'■^1 >-^2 A 


— A(mi+m2) 


i——t m2 ——t2 
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hence it is sufficient to calculate ^^(—2). We will show by induction on n that indepen¬ 

dent of {n,k), or equivalently that ^2(~2) is independent of {ii,£ 2 )- Since we have established the 
case n = 0, Lemma [5]8] then follows. 

In order to perform the induction step we consider the recursion of Theorem 14.61 Using = 

we see that </j(1) = + q~^). Taking the trace of Theorem 14.61 at A^ = 1 

and using Proposition 15.2f iil we find, after taking traces, 

(6.4) '^{q + q + 1) = (4ll/2,l/2 + 4l_i/2,-l/2 + ^- 1 / 2 , 1/2 + £^ 1 / 2 ,- 1 / 2 ){‘^£ + !)■ 

Next we evaluate Theorem 14.61 at A~^ and we take traces, so, using ip(A~^) = 

2(9 + 9 ^)//i,^2(“2) = X/ 

iJ=±l/2 

which we rewrite, assuming ~ independent of {£i,£ 2 ) for li -\- £2 < £ n, as 

/ 4 +i,^ 2 +i (-2) = (^ 2 *''^ ~ M/ 2 - 1/2 - -4-1/2, 1/2 - -4-1/2,- 1 / 2 ^ = Ft, 

by (16.41) for the last equality. So the statement also follows for £1 -|- £2 = f -I- n -I- 1, and the lemma 
follows. □ 


7. g-DiFFERENCE OPERATORS FOR THE MATRIX-VALUED POLYNOMIALS 

We continue the study of g-difference operator for the matrix-valued orthogonal polynomials started 
in Corollarv l5.13l In particular, we show that the assumption on the invertibility of $0 follows from the 
LDU-decomposition in Theorem 14.151 Then we make the coefficients in the matrix-valued second order 
(/-difference operator of Corollary 15.131 explicit in Section 17.11 Comparing to the scalar-valued Askey- 
Wilson (/-difference operators, see e.g. [2], [16], [20], [26], we view the (/-difference operator as a matrix¬ 
valued analogue of the Askey-Wilson (/-difference operator. Next, having the matrix-valued orthogonal 
polynomials as eigenfunctions to the matrix-valued Askey-Wilson g-difference operator, we use this to 
obtain an explicit expression of the matrix entries of the matrix-valued orthogonal polynomials in terms of 
scalar-valued orthogonal polynomials from the g-Askey scheme by decoupling of the g-difference operator 
using the matrix-valued polynomial L in the LDU-decomposition of Theorem l4.15l From this expression 
we can obtain an explicit expression for the coefficients in the matrix-valued three-term recurrence 
relation for the matrix-valued orthogonal polynomials, hence proving Theorem 14.111 and Corollary 14.121 

7.1. Explicit expressions of the g-difference operators. In order to make make Corollary 15.131 
explicit, we need to study the invertibility of the matrix $q(A'''). For this we first use Theorem 15.51 and 
its Corollarv l5.61 and in particular (16.11) . Because of ProDOsition l5.2l iil. this is only a single sum. In the 
notation of (I5.15L we find 

i 

W(iA)fc.p(A^) = ^(4.^)„,,(A^-i)($g)„,^(A->-i) = ((ci.^(A-^-i)) 4>^(A^-1)) 

n=-t 

Taking the determinant gives, recalling /i(z) = ^{z + z~^), 

det (dg(A-^-i)) det = det(lU(/r(g^))) = det(r(/i(g^))) 

2t 21 

= n r(/r(g^))fe,fe = ^4cfe(^)(9'+'^9"■'^9')fc 
fc =0 fe =0 

using the LDU-decomposition for the weight of Theorem l4.15l and (14.121) . The right hand side is non-zero 
for A G Z unless 1 < |A| < 2£. So <i>Q(A'’') is invertible for A > 2f or A < —2£ — 1 or A = —1, i.e. for an 
infinite number of A G Z and it is meaningful to consider Corollary 15. 131 

Proof of Theorem \4-l^ Corollary 15.131 gives a second-order g-difference equation for the matrix-valued 
orthogonal polynomials for an infinite set of A. So it suffices to check that for z = 1, 2 the matrix-valued 
functions Ali, Ni of Corollary 15. 1 31 coincide with Mi, Nt, where Afi{z) = Mi(z~^), of Theorem 14.131 or 

(7.1) $^o(A^-i)7W,(g^) = M,(g^-i)4.^(A^), ip^M(9^) = 

where Mi, Ni as in (I5.16p . see Proposition l5.11[ and Mi, Ni as in Theorem l4.13[ where M(-2^) = Mi{z~^). 
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By (I5.18P we need 


KQ{i) = Mi{z) + Mi{z) = Mi{z) + Mi{z ^), 


(7.2) 

which is an easy check using the explicit expressions of Theorem 14.131 Now (17.2|) and (15.161) for n = 0 
show that any equation of dm implies the other equation of dZH)- Indeed, assuming the second equation 
of (17. 1|) holds, then 

= $g(7l^)Ao(z) = + N,{q^)^iA^-^) 

= M,iq^)^{A^+^) + ^{A^)M,{q^+^) 


implying the first equation of dUD- 

By ProDOsition l5.2l iil and (I5.15|) the matrix entries of <I>q(A'’') are Laurent series in q^. Setting z = q^, 
we see that in order to verify dm entry-wise, we need to check equalities for Laurent series in z. 

We first consider the second equality of dm for f = 1. In this case the matrices iVi and Ni are 
band-limited. Hence, the (to, n)-th entry of both sides of dm involves either two or one terms, so we 
need to check 


(7.3) 


^ A/*! (z)n-l,n + ^^(o,n) )m-e,m-£ _ 


z—q-^ 




The proof of (17.31) involves the explicit expression of the spherical functions in terms of Clebsch-Gordan 
coefficients using Proposition 15. 2l iil. It is given in Appendix IB. 31 

The statements for the second g-difFerence equation with i = 2 follows from the symmetries of Propo¬ 
sition 15.111 and Lemma 15.141 □ 


The explicit expressions have been obtained initially by computer algebra, and then later the proof 
as presented here and in Appendix IB. 31 has been obtained. 


7.2. Explicit expressions for the matrix entries of the matrix-valued orthogonal polynomials. 

Having established the (/-difference equations for the matrix-valued orthogonal polynomials of Theorem 
14.131 and having the diagonal part of the LDU-decomposition of the weight in terms of weight functions 
for the continuous g-ultraspherical polynomials in Theorem l4.15l it is natural to look at the (/-difference 
operators conjugated by the polynomial function L*. It turns out that this completely decouples one of 
the second order (/-difference operators of Theorem 14.131 This gives the opportunity to link the matrix 
entries of the matrix-valued orthogonal polynomials to continuous g-ultraspherical polynomials. In order 
to determine the coefhcients we use the other (/-difference operator and the orthogonality relations. 
Having such an explicit expression we can determine the three-term recurrence relation for the monic 
matrix-valued orthogonal polynomials straightforwardly, and hence also for the matrix-valued orthogonal 
polynomials P„, since we already have determined the leading coefficient in Corollary 14.91 

The first step is to conjugate the second order (/-difference operator Di of Theorem 14.131 with the 
matrix of the LDU-decomposition of Theorem 14.151 into a diagonal (/-difference operator. This con¬ 
jugation is inspired by the result of [Ml Theorem 6.1]. This conjugation takes in a three-diagonal 
(/-difference operator. For any n G N, let TZn{x) = L*{x)Qn{x), where Qn{x) = Pn(a;)(lc(P„)) ^ denote 
the corresponding monic polynomial. Note that we have determined the leading coefficient lc(P„) in 
Corollary 14.91 Then (TZn)n>o forms a family of matrix-valued polynomials, but note that the degree 
of TZn is larger than n, and that the leading coefficient of 77.„ is singular. Note that 77.„ satisfy the 
orthogonality relations 


(7.4) 


r*l _ /*1 

(TZn{x))*T{x)'R-m{x) \/l- x^dx = J {Qnix))*W{x)Qmix) V^ - x'^dx 


= 5m,n\{\c{Pmr) 'G„(1c(P™)) \ 


Theorem 7.1. The polynomials (7^n)n>o ctre eigenfunctions of the q-dijference operators 
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with eigenvalues A„(i), where 


It. 




i=0 


(l-g2)2 (l-z2) 


2^ 


n _ „4^-2i+2'| 
i-2^+1 U g ) z 


K.,(A = 


i=l 

2^ 


1 


(1 + g4^+2) (1 - g2*+2^2) 


E 0„®-2^+1 ~ _ ' \~ ~L ~ / p 

^ (l-g2)2 (l + g2,)(l+g2»+2) (1-^2) 


2^-1 


V- *-2m 1 (l_g4£+2.+4)(i_^2.+2)2 _ ^2*+2^2) _ ^2.+4^2) ^ 

“ 2^ 9 T:--I-n^TT^-T-71TTT^-- o.-,o^o -7^-^7-A 


i=0 


(1 — g2)2 _ g4j+6)^2. — g4*+2)(J^ _|_ g2i+2^2 


Z(l-Z2) 


z,i+l- 


Proof. We start by observing that the monic matrix-valued orthogonal polynomials Qn are eigenfunctions 
of the second order g-difference operators Di of Theorem l4.13l for the eigenvalue lc(P„)A„(i)lc(P„)“^ = 
An{i), since the matrices are diagonal and thus commute. By conjugation we find that satisfy 

ICi{z)iin{qz) + ICi{z~^)iin{^) ='R.n{x)An{i), ICi{z) = L\z)Mi{z) (L*(qz)) ^ 

using the notation P‘(z) = L*{g,{z)), etc., with x = ^{z) = ^{z + z~^) as before. It remains to calculate 
JCi (z) explicitly. We show in Appendix IB. 41 that the expressions for fCi are correct by verifying 

(7.5) }Ci{z)L*{qz) = L*{z)Mi{z), i = 1,2. □ 


Lemma 7.2. For n G N and 0 < i,j < 2£ we have 

'R-n{x)ij = /3„(i,j) C'„+j_i(a:;g^*+^|g^), 

where Cn{x] P\q) are the continuous q-ultraspherical polynomials (14.101) and Pn{i,j) is a constant de¬ 
pending on i,j and n. 


Proof. Evaluate 'DplZnix) = 77.„(x)A„(l) in entry {i,j). Since Vi is decoupled, we get a g-difference 
equation for the polynomial (77.„). ., which, after simplifying, is 


(1 _ 52 *+ 2^2 
(1-Z2) 


-iin{qz)i 


(1 — g^*+^z 


P-n(g 2)ij=g+*(g 




+ g^'+"+i)P„(z),,. 


All polynomial solutions of this g-difference are given by a multiple of the Askey-Wilson polynomials, see 
pm §16.3], [ini §7.5], pni §14.1], p„-i-j-i(a;;g*'*'^, —g*+^, g^/^, —g^/^jg). Apply the quadratic transforma¬ 
tion, see [U (4.20)], to see that the polynomial solutions are Pn+j-i{x-, g*+^, —g*+^, g*+^, —g*'''^|g^). These 
polynomials are multiples of continuous g-ultraspherical polynomials Cn+j-i{x', g^*+^ |g^), [20l §13.2], [161 
§7.4-5], pni §14.10.1]. Hence, 'Rn{x)ij is a multiple of C'„+j_i(a:; g^*+^|g^). □ 


Our next objective is to determine the coefficients fin(i,j) of Lemma [7.21 Having exploited that the 
matrix-valued polynomials 77.„ are eigenfunctions for the decoupled operator (Di of Theorem 17.11 we can 
use Lemma O in (17.41) to calculate the (i,g)-th coefficient of (17.41) : 

t -7 a\ -^Pn{k,i)Pm{k,j)ck{£) f C„+j-fc(a:;g^^+^jg^jCrn+j-fc(a:;g^^+^|g^) ^ dx 

( 7 - 6 ) J-i Vl-x^ 

(G 

m )^,^(lc(P„))-^ 

using that lc(Pm) and the squared norm matrix Gm are diagonal matrices, see Corollarv l4.9l and Theorem 
EM The integral in (17.61) can be evaluated by (14.111) . In particular, the case m + j = n + i of (17.61) gives 
the explicit orthogonality relations 

(9 ,9 )fc (9,gOm-Hj-fe (1 - g2m+2j-H2) 

„ (i-g4^+2)2 (g^g4'+^g2)^ 




_ g2?7T.-[-22-1-2 ^ ^ ^ — 2 2-j-2 771-[-2 


) (g^*+^g 


,4^-22-H2. ^2^12 ■ 
’ y )m 
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Theorem 7.3. We have 


(„2 M+i. .,2\ «- 2 i- 2 rt. „ 2 'i (^2. „2\ 

-U —{ lUo-rt W ig 1 g jri Ig ig _ Jz Ig ; g M+j-z ,-(2i+l)+2i(2^+7t+l)--i^ 


Proof. From Theorem 17.11 we have 

(7.8) lC 2 {z)n^{qz) + lC 2 {z-^)nn{q-^ z) = n^{z)Kn{2). 

Evaluate (TTSl) in entry (i,j) and use Lemma [7.21 to find a three term recurrence relation in i of /3„(z, j), 


(9 


-j-n-l 


q- 


j+n+l 


(q -1 - 9)2 


Pn{t,j)Cn+,-^{z;q^^+^\q^) 


= j3n(i + 1 , j)(^/C2(z)iy+iC„+j_i_i(9z;9^*+"^|9^) + K.2{,z~^)^^i+lCn+j-i-l{q~^z]q^'^^\^ 

+lin(i, j){r 2 {z)uCn+i-i-i{qz-,q^''^‘^\q^) + K,2{z~^)uCn+3-t{q~^ Z\q^'^^‘^\q^)^ 

+/3„(z - 1, j), (^/C2(z)i,*_i(7„+j_i+i(9z; 9^*|9^) + K.2{z~'^)i^i-iCn+j-i+i{q~^z', 9^*|g^)) 

Multiply by (1 — z^)(l — 9 ^)^ and evaluate the Laurent expansion at the leading coefficient in z of degree 
n + 9 — z + 3. The leading coefficient in z of the continuous 9 -ultraspherical polynomial Cn{az] l3\q) is 
(g’g)” After a straightforward computation this leads to the three term recurrence relation 

(l + g"'-"^'+2"+2)/3„(z,9) = 

2.+3 (1 - 94^+2»+4)(1 - g 2 »+ 2 ) 2 (l - g2n+2z+2,-+4) 

^ (l-9"+*+J+2)(l_g4»+6)(l_^4*+2)(l_^2*+2)2/^"V*+ ,J) 

2,+2 (1 + 94 ^+ 2)(1 - g 2 «+ 2 j+ 2 ) 


- 29 " 




(l + 92*)(l + 92i+2) 

+9^*+^(l - 9"+*+-^+i)(l - 9'^^“^*+2)(1 - 92 "+ 2 j- 2 *+ 2 )^^(j _ 1 ^ 


This recursion relation can be rewritten as the three-term recurrence relation for the 9 -Racah polyno¬ 
mials after rescaling, see [TCI Ch. 7], [TCI §15.6] [TCI §14.2]. We identify {a,/3,j,S) as in (14.131) with 
( 1 , 1 , 9 “ 2 «- 2 i- 2 ^ ^- 4 £- 2 ^ l^g^gg ^ 2 _ gives 


Pnii,j) 


7n(j)(-l) 


, (9 ‘^^q ^"; 9 ^)^ {q‘^-,q^)u+j-i 2n+2i(2e+n+l)-i^ 

( 92 , 94 ^+ 4 ; 9 ^), {q^^+^;q^)n+J-^^ 

xR,{n{j); 1 , 1 , 9 “^^"^; 9 ^) 


for some constant 7 ^( 4 ) independent of z. Plugging this expression in (17.711 for i = j gives l7n(9)l^ by 
comparing with the explicit orthogonality relations for the 9 -Racah polynomials, see [M Ch. 7], [TCI 
§15.6] [TC] §14.2]. 

For j > i we have Tiij{x) = Lj_i(x)(x"Id -I- l.o.t), and since the explicit expression of Ljy shows that 
the leading coefficient (of degree j — i) is positive, we see that the leading coefficient (of degree n + j — i) 
of TZij in case j > z is positive. Since 'Jnij) is independent of z, we take z = 0, which shows that 7 n(j) 
is positive. □ 


Proof of Theorem \4-17\ Using Theorem 17.31 with the explicit inverse of L{x) as given in Theorem 14.151 
gives an explicit expression for the matrix entries of Qn{x) = (L(x)“^)*77.„(x). Then we obtain the 
matrix entries of Pn{x) = Qn{x)\c{Pn) from this expression and Corollary 14.91 stating that the leading 
coefficient is a diagonal matrix. □ 


7.3. Three-term recursion relation. The general matrix-valued orthogonal polynomials satisfy a 
three-term recurrence relation, see Section |2| Moreover, Theorem SJ shows that the three-term recur¬ 
rence relation can in principle be obtained from the tensor-product decomposition. However, in that 
case we obtain the coefficients of the matrices in the three-term recurrence relation in terms of sums of 
squares of Clebsch-Gordan coefficients, and this leads to a cumbersome result. In order to obtain an 
explicit expression for the three-term recurrence relation as in Theorem 14.111 and Corollary 14.121 we use 
the explicit expression obtained in Theorem 14.171 and Lemma FfTC which is [TC] Lemma 5.1]. Lemma [7TC 
is only used to determine X„. 
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Lemma 7.4. Let (Qn)n>o be a sequence of monic (matrix-valued) orthogonal polynomials and write 
Qn{x) = , where G Matjv(C). The sequence {Qn)n>o satisfies the three-term recurrence 

relation 

xQnix) — Q-\- CfjlixfjXji -j- 
where Y-i = 0, Qei(x) = I and 

Y _ _ /n«+l V _ D”' _ _ G)" Y 

— 'Yn—1 Vrt ) ^ n — ^^7^—2 ^n — 1 Wn — l^n- 

So we start by calculating the one-but-leading term in the monic matrix-valued orthogonal polynomi¬ 
als. 


Lemma 7.5. For the monic matrix-valued orthogonal polynomials {Qn)n>o we have 


21.-1 


^q (l-g2n)(i_g2,-+2) (1 _ _ ^4£-2,-+2) 

^ 2(l-g2)(l_g2n-h2j+2) 3 J + 1 ^ 3 (1 - g2) (1 _ g«-2j-S2n-h2) JJ-l’ 

Proof. Since {TZn{x))ij = and Qn(x) = (L*(a;)) ^Tln{x) by Theorem l7.3l see 

Section [El we have 

(7.9) Qr^ixfi,, ^ £^( 2 fc+l)(fc-i) (g ;g h(q ;q y 

W I 9 ) 2 k[q ; <? ji 

and this expression shows that deg((5„(a;)ij ) = n -I- j — L So in case j — i < 0, we can only have a 
contribution to Qn_i in case i — j = 0 or z — j = 1. The first case does not give a contribution, since 
(17.91) is even or odd for n -I-j — z even or odd. So we only have to calculate the leading coefficient in (17.91) 
for z — j = 1. With the explicit value of /3„(z, j) as in Section |E1 or Theorem 14.171 and Corollary 14. 19l we 
see that Q"_i in case z — j = 1 gives the required expression for (Qn-i)^ 

On the other hand, by Proposition 14.101 and since J\c{Pn)J = ^c{Pn) by Corollary 14.91 it follows that 
JQn{x)J = Qn{x). This gives (Qnix)).^ = (Qnix)) .^g_. 2 t.-p t^ 2 at the case j — z > 0 can be reduced 
to the previous case, and we get {Ql-i) ^ ° 

Proof of Theorem \4.11\ The explicit expression for follows from Lemma 17.41 and Lemma 17.51 
By Theorem 14.81 and Corollary 14.91 we have the orthogonality relations 


n 1 

J {Qm{x))*W{x)Qn{x)\/^ - x'^dx = Sm,n{^c{Pm)*) ^Gm(lc(P„^)) ^ = Sm,nGm{^c{Pm)) ^ 


since the matrices involved are diagonal and self-adjoint, hence pairwise commute. By the discussion in 
Section H we have 

Y„ = G-^(lc(P„_l))'G„(lc(P„))■^ 

and a straightforward calculation gives the required expression. □ 


Appendix A. Branching rules and Clebsch-Gordan coefficients 


Proof of Theorem \4.S\ The Clebsch-Gordan decomposition for Uq{su{ 2 )) and the corresponding inter¬ 
twiner involving Clebsch-Gordan coefficients is well-known, see e.g. [21] §3.4]. With the convention of 
the standard orthonormal bases as in Section [3] the Clebsch-Gordan coefficients give the unitary inter¬ 
twiner 

= 0 iLifi- 0 

fc|A-A| fc|A-A| 


(A.1) 




ha V V pA.A.t Ji ^ J2 

^ / J / J ^ni,n 2 ,n ^ni ^ ^n 2 ’ 

ni =—n 2 =—^2 

-^l +-^2 

^ t\X), VAgW,(su(2)). 


Identifying the generators ^ x x, F) of the Hopf algebra as in [3TJ §§3.1.1-2] with the generators 
(k-G\lfil\f^e) as in Section 131 we see that the Hopf algebra structures are the same. Moreover, in this 
case the representations fi' as in Section [S] correspond precisely with the representations T\i of [HI §3.2.1] 
including the choice of orthonormal basis. This gives = Gq(£i, € 2 ,n-i, n- 2 , n-), where the right 
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hand side is the notation of the Clebsch-Gordan coefficient as in [5T1 §3.4.2, (41)]. The Clebsch-Gordan 
coefficients are explicitly known in terms of terminating basic hypergeometric orthogonal polynomials, 
the so-called g-Hahn polynomials m , see m §3.4]. 

In order to obtain Theorem l4.3l from we use Remark l4.2l iil. With the ^-algebra isomorphism 'k 

as in Remark 14.21 we have = J^t^{X)J^ for all X € Wq(su(2)), where the intertwiner J^: —?> 

Cp !->■ elp is a unitary involution. Note that the representations of Uq{su{ 2 )) in (13.3|) and of B 
in (14.21) are not related via the map of Remark I4.2l iil. but they are related by the same operator J^. 
Theorem now follows after setting G Id) o 7 ^^ o J^, so that in particular = 


^£1 ,£2 

ni ,712 ,—n 


-ni,n 2 , -n). 


□ 


The Clebsch-Gordan coefficients satisfy several symmetry relations, and we require, see [HJ §3.4.4(70)], 

/A o'! 

^ni,n 2 ,n ^n 2 ,ni,—n' 

We need explicit expressions of the Clebsch-Gordan coefficients for the case -|- ^2 = which follow 
from the explicit expressions in EB §3.4.2, p. 80]. For fixed t S ^N, let —£< m < £, and we consider 
the case £i = {£ + m)/2 and £2 = {£ — to)/ 2. The Clebsch-Gordan coefficients in this case are given by 


(A.3) 






£-\-m 

g2 

l—m 



21 

£-k 



assuming i — j = k and i £ {—^i£ + rn), ^{£ -|- m) -I- 1, • • • , ^{£ + m)}, j £ {—^{£ — rn), ^{£ — m) + 
!,■■■ ,4(£-m)}, k £ {-£,-£+ 1 , • • • ,£}. 

Observe that the unitarity gives 


(A.4) 


X — Ip^ P^\ — fill,£2.^ fill,£2/ i ll pl2 pll (^p£‘ 2 .\ 

ni,mi=— n2,rfi2 ——l2 

li I 2 

_ \ ^ \ ^ fill,£2,1 fill,£2,1 

/ > / > ^ni,n2,m^ni,n2,n 

ni——£i n 2 ——£2 


using the fact that the Clebsch-Gordan coefficients are real, see EH §3.2.4]. 

We also need some of the simplest cases of Clebsch-Gordan coefficients, see ED §3.2.4, p. 75], 


(A.5) 


C. 


1 / 2 , 1 / 2,0 


-1 


c 


1/2,1/2,0 


t/ 2 , 1 / 2,0 yrr^’ -- 1 / 2 .- 1 / 2.0 
1 /2 1/2 0 

and the other Clebsch-Gordan coefficients q ’ being zero. Another case required is a generalisation 

of a case of (ESI) , see ED § 3 . 2 . 4 , p. 80 ] 


(A. 6 ) 

and more generally 


_ 21 
'~'-e,-e,o ~ 9 1 


' l-g 2 

1 - g4^+2 


(A.7) 


C‘ 


'll ,£2 ,l 


— ll, —12 ,l2~ll ^ 


_ ^£i-\-l2~l 




2^2 


(9^; 9^)4+A+f-ei(9^; q'^)ei+e2-i 




Appendix B. Proofs involving only basic hypergeometric series 

In Appendix|B]we collect the proofs of various intermediate results only involving basic hypergeometric 
series. For these proofs we use the results of Gasper and Rahman m- In particular we follow the standard 
notation of [16], and recall 0 < g < 1 . 

B.l. Proofs of lemmas for Theorem 14.81 Here we present the details of the proof of Lemma [6.21 
which is used in order to prove the explicit expression of the matrix entries of the weight in terms of 
Chebyshev polynomials in Theorem 14.81 We start with two intermediate results needed in the proof. 
Lemma IBTI can be viewed as a g-analogue of Sheppard’s result [21 Cor. 3.3.4]. 
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Lemma B.l. For b, c,d,e € we have 


Y-/, fe N / fe \ (9 '',b,c]q)k (de , 

y{dq ]q)n-k[eq ;q)n-k -7-^- -j-g 


/c=0 


- 2 (S) 


)n — k " 


(<? ",c,bcq^ '^/{de);q)k ^ 


( 9 ; 9 )^ 


-9 • 


Proof. Applying first [1^1 (III.13)] and next [THl (III.12)] gives 

fq-'^,b,c de {elc\q)n fq-'^,c,d/b 

^ „ (d/c, elc\q)n /g“”,c,6cgi"”/de_ g 
^ id,e;q)n jd,cq^-^jb 

Multiplying with (d, e; g)„ and expanding gives 

V-/J fc k ^ {q~'^,b,c]q)k fde {dlc,efc;q)n {q-^,c,bcq^-'^/de-,q)k (q\^ 

y,(d, ,e,(fcd j =■= L 






fc =0 

(g-"',c,6cgi-”/(de);g)fc ^ 


{q]q)k 


(I)' 


)n — k " 


{q;q)k 


-q ■ 


□ 


Lemma IB . 2 1 Efives a simple g-analogue of a Taylor expansion. 


Lemma B.2 (a g-Taylor formula). If we have a polynomial B{q ^] (l)t, then 


KB{q-^) 


1\ 


— Ari 


M =0 


where Aq is the q-shift operator AqB{q “) = q^{B{q ”) — B{q “ ^)). 

Proof. Define ft{x) = {x\q)t, then 

^Jt{q-^) = q^iftiq-^) - Mq-^-^)) = q{l - q*){q-^; q)t-i = g(l - g‘)/*-i(g""^). 

Repeated application of Aq on ft{q~^) then gives 

(B.l) A^/t(g-^) = g"(g‘;g-i)„/t_„(g-“) =g-+‘"-(S)(-l)"(g-‘;g)„/*_„(g-^). 

Putting M = 0 in ()B.lll the expression is zero if n > t and n < t. Therefore 

This gives the result. 


□ 


Proposition B.3. Let £ S and k,p gN such that k,p<2£, k — p<0 and k + p < 2£. Take s G N 
such that s < p and define 


esik,p) =q 


p—s 

_ ^2ks-\-k-\-p—2i E g-2i(s+l) 


2=0 


'k' 


p 

i 

q2 

i + s_ 


k—s 

n=0 


Ai-k 


2i-p 

n+s 


n 


2e 

2 +n+s 


1 ^ 
- 


Then we have 


£(, N _ 2v(2i+i)-p^-2i+k (1 - {q^-,q%i-piq^;q^)p 

« (l-g2fc+2) (g2;g2)2, 




p-T 


T=0 


(g^^-^^;g^)p-T (g4^-^^+"; g^), 

(g2fc+4;g2)p_5, {q^\q^)T 
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Proof. The proof proceeds along the lines of the proof of [531 Proposition A.l], and we only give a sketch. 
We start by reversing the inner summation, using M = 2i — k — s — m in el{k,p) to get 


(B.2) 


e^(fc p) = g2fcs+fc+p—2^+2(fc+p—s+l)(2f—fe —s) ’'y ^2i(2fc+s —1) 

) 

[ 


2=0 


'k 


p 

i 

g 2 

.i + s_ 


J 


k—s 


2l-k 

2M{2i-k-p+s-l) ^ iq^ 


21-p 

2i — k — s — M 


M=0 

We rewrite the inner summation over M 

k—s 


2t 


1 ^ 

- 9^ 




2i-k 2i-p 

, .Mo 2£—k — s — M 

2M{2i-k-p-\-s-l) L J I 


M=0 


(B.3) 


2i 

2i—p—M-\-i 


1 ^ 
- 


(g^; q^)2i-k{q^\q^)2^-p{q^]q^)2i-k+^(.q^]q^)k-^ 


{q^;q^)Uq^;q^)2i-k-s 

^ 2M(k+s) g^)M p. -2 Mn 


M=0 


{q^\q^)M{q~ 


*;g 2 ) 


where B{q~‘^^) = is a polynomial in g"^'^ of 

degreep—s that depends on I, k,p, i and s. The polynomial B{q~‘^^) has an expansion in (—l)*(g“^'“; q‘^)t 
such that B{q~'^^) = X]t=o At(—l)*(g“^'^; g^)t. By Lemma lR2l the coefficients At are obtained by re¬ 
peated application of the g-shift operator; 

ASB(g- 2 ^) 

(B.4) -— = A*. 


g2(‘t')(g-2t;g2), 


M=0 


We substitute B{q~‘^^) = At(—l)*(g“^'^; g^)t into (IB.3|) and we interchange the summations over 
M and t. Then the sum over M can be rewritten as a summable 2 <Pi- Using the reversed g-Chu- 
Vandermonde summation [131 (II-7)] we can rewrite (IB.3I) as 

., 2 .„ 2 \ /„ 2 .„ 2 ^ l„2.„2\ l„2.„2\ 


{q ]q ) 2 e-k{q ]q ) 2 e-p {q ; q%i-k+i{q ■, q )k-i 


(B.5) 


(g2;g2)2 


2e 


{q^\q^)2i-k-a 

^ _2t(fc+s+t+l)+2(‘) . (g2fc-4^+2«^ g2fc-2»-e2. ^2)^(g-4^-2. 

‘ (g-«+ 2 fe- 2 ,.^ 2 )^^^^ 

Substituting (M and (IBAI) into (IB. 21) and interchanging the sums over i and t, we obtain 
Z(h A 2ka+k+p-2e+2(2e-k-s)(k+p-s+l) (g^ g")fc(g^ g")p (g" ( g" )2^-fc (g^ g")2^-p 


eaik,p) = g 


(g2; g2)2^(g2; g2),,_,_,(q2; ^2)^(g2; 

/^^-2*(fc+«+*+l)-2(‘)(^2fc-4.+2s.^2)^(g-4.-2.^2)^^_^_^_^(^2.+2.^2)^ 


p—s 


t =0 


Q 


g 2 (‘t')(g- 2 *;g 2 ), 


q 


2M{s—p) ^^2k — 2p-\-2s-\-2M-\-2. ^2 


\q%-a 


M=0 


p—s 


X y^^q2t{-2i+p-s-t-l) 


i=0 


(„2s-2p „-2fc „it-2k-2M+2. ^,2\ 

\q _ ^q A _; q )i 

(g^, g^)i 


The inner sum over * is a 3 P 2 -series and after some rewriting, it can be transformed using Lemma |B. II 
The sum over i becomes 


p—s 

E' 

i=0 


2i{-2l-\-p-s-t-l) ^^-2k-2M-\-2i ^-2k-2t+2i. ^2 




) Q. )p—s—i ■ 


/„4^-2fc-2M+2 ^-2fe „2s-2p.„2\ 

Vi o y iH : H )i 


iq^;q^Z 


= q 


p—s 

Er 

2=0 


... —2p+2s „— 2k ^4i£-\-2t —2p+2s+4. „2'\ 

-2A:(p-s) ^22(2p-2s+fc-t-2^-2)-4(^)^^-2i ^-2M.^2', W i Q ^ Q iH )i 


j q )p—s—i 


{qZq^Z 
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We observe that the evaluation of the t-order g-difference operator yields only one non-zero term in the 
sum over i, namely the one corresponding to i = p — s — t. After simplifications we have 

_ J4l-2v+3)s+s^+k-2i-v (1 - (g^; q^)2l-p{q^]q% 

® (l-g2fc+2) 


p—s 


E(-i) 


, („2k-U 2\ ( Al-2p+2s+2t+4 2\ 

s+t t^+t{4e-2p+2s+3)^ 5? ,q )p-s-t 


t=0 


{q'^>=+'^;q'^)s+t{q'^-,q'^)p-s-t 
Reversing the order of summation using T = p — s — t proves the proposition. 

Proof of Lemma \6.2\ Recall from Proposition 16.11 that we have 

\(k+p) 

^ di{k,p)q^^\ 

s=-4(k+p) 

Therefore the coefficients of dl{k,p) and a^{k,p) are related by 

p-(s+|(p-fc)) 

diik,p) = (A(k,p). 

t=o 

Suppose that p — k <0, k+p <21 and s > \{k — p). Using Proposition 16. II we find the expression 


□ 




_ 2ks-\-k-\-p—2t 


p—s 

E' 

2=0 


2i(s-|-l) 

'k 


p 


i 


i -(- s_ 


2i-k-s 

E 

9^ n=0 

it 


q 


2n( — 22 +fc+p—s+1) Y 


2£-k 

n+s 

g2 

2e-p 

n 



2£ 

2 


By taking ei{k,p) = . i,. „ Proposition IB. 31 shows that , i.. Ak,p) is equal to 

^2p{2i+i)-p^-2t+k (1 - {q^]q^)2t-p{.q^-,q^)p 

(B. 6 ) 


(I_g 2 fc-H 2 ) {q^]q^)p 

p-(4t+3)t, . ^k-t q^)p-t (^4^-^*+-^ ; q% 

^ ^ (g2fc+4;q2)p_^ (^2.q2 


p—s 


t =0 


Comparing (IB.61) with the explicit expression of a^{k,p) yields the statement of the lemma. 


□ 


B.2. Proofs of lemmas for Theorem l4.15L Here we present the proof of the technical lemmas needed 
in the proof of Theorem 14. 151 in which the LDU-decomposition of the weight matrix is presented. 

For completeness we start by recalling the linearisation and connection relations for the continuous 
g-ultraspherical polynomials, see [H §10.11], [101 §13.3], (TOJ (7.6.14), (8.5.1)]. The connection coefficient 
formula is 


Ln/2J 


(B.7) 


Cnix-,-f\q) = Y 


1 - Pq 


n — 2k 


A ;=0 


1-/3 


jfc {i/P;q)k{i;q)u-k 

{q;q)kiPq; q)n-k 


C„-2k{x;P\q), 


and the linearisation formula is 


(B. 8 ) 


C„(x;/31g)C™(x;/31g) = Y ^ iq,qU+n- 2 k{P; qU-kiP', q)n-k 

1 P \P ^ Q)m-\-n—2k{Q^ Q)m — k\Q'i Q)n—k 


fe =0 


iP;q)kiP'^;q) 


m-\-n — k 


Cm-\-n—2k(^^i P\q')- 


{q\q)k{qP]q)m+n-k 

The proof of Lemma 101^ follows the fines of the proof of [HI Lem. 2.7] closely. 

Proof of Lemma \6.4\ In order to evaluate the integral of Lemma l6.41 we observe that the weight function 
in the integral is the weight function (14.1111 for the continuous g-ultraspherical (in base q^) with /3 = . 

Rewrite the product of the two continuous g-ultraspherical (in base g^) with /3 = g 2^+2 ^ sum over i 

of Cm-i-n- 2 fc- 2 i(a:; g^^'''^lg^) using (|B. 8 p . Since the Chebyshev polynomials can be viewed as continuous 
g-ultraspherical polynomials, Um+n-2t(x) = Cm+n-2t{x] g^Jg^), we can use (|B.7|) to write the Chebyshev 
polynomial as a sum of continuous g-ultraspherical polynomials with P = Plugging in the two 
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summations, and next using the orthogonality relations (14.111) shows that we can evaluate the integral 
of Lemma inU] as a single sum; 


/^2fc+2 „2fc+4. „2'\ 
( 2 fc+ 2 )(fc-t) W jQ )c 


(q,2^g4fe+4.q2)^ 


(B.9) 


tAm — k 

E 

r—OVt—k 

X 


_ ^2r»+2n-2fc+2-4.) (g2fc+2.g2)^ ( 9 ^^+^ ; 

(1 _ q 2 m+ 2 n- 2 fe+ 2 - 2 r) (g 2 .q 2 )^ {q^; q'^)^_k-r {q^; q‘^)n-k-r 


/'„4fc+4. _2'i t„—2k.„2\ („2..,2\ 

\Q )m+n—2k—r [Q j <? J/c+r—t V9 )m+n—t—k—r {2k+2)r 

We consider two cases; k > t and k < t. li k > t, note that 

^ g2m+2n—2fc+2—4r 


= 9 


m—n+1 m—n+1 „2fe—2m—2n—2. „2\ 

^_2r W 5 y 1 y iH )r 

(nk—m—n — 1 nk—m—n — 1 ri‘^k — 2m—2n.r,2\ ’ 

W 5 y 5 y 'iH )r 


_ q2m-\-2n—2k-\-2—2r 

so that we can rewrite (El as a terminating very-well-poised g^/^T-series. Explicitly, 

(g2fc-^2^g2fc+4;g2)^ (g2fc+2;g2)^_^ (g2fc-^2, ^2 (g^; 

(^q2^q4k+4.q2'^^ (q2;g2)^_^ (q2;g2)^_^ ((j2fc+2; (q2;g2)^_^ ((j2fc+4;g2)^_^^_^_^ 


X g(2fe+2)(fe-i) 


2Ai—2m—2n —2 2fe —2m 2k—2n 2t — 2m—2n—2 2k-\-2 2 —2k\ 

?y ^y ^y ^y ^y •> Q ■i Q )• 


Here we use the notation of [m § 2 . 1 ], Using Watson’s transformation formula [THl (III. 18)] and recalling 
the definition (I4.13|) of the g-Racah polynomials we can rewrite the sll? as a balanced 4 'Ci 73 ; 


2 fc- 2 t -|-2 q-2m-2n-2. q2 

/„ 2 fc— 2 m— 2 n „— 2 i, „ 2 '\ 

_ W _ ? 9 ; g jt 

^g 2 fc— 2^+2 q—2m—2n—2.q2'j^ 


-2k „ 2 i- 2 m- 2 n -2 ^ 2 A :+2 ^- 2 t 

t /y?y ^y ^y 22 

4^3 ( q-2m^q-2n^q2 


(^2k—2m—2n „— 2 i. „ 2 \ , „ ^ „ 


Rfc(Mi);i,i,9-“,9-""-";9"). 


Simplifying the (/-shifted factorials gives the required expression. 

For the second case k < t, we proceed similarly. After applying Watson’s transformation we also 
employ Sears’ transformation for a terminating balanced 4(/?3 series |161 (III.15)] in order to recognise 
the expression for the (/-Racah polynomial. □ 

We leave the proof of the generalisation of Lemma 16.41 in Remark 16.51 along the same lines to the 
reader. 

In the proof of the LDU-decomposition of Theorem 14.151 we have also used Lemma 16.61 whose proof 
is presented next. 

Proof of Lemma 1 6‘. 61 We Hrst write the g-Racah polynomial in the left hand side of Lemma 16.61 as a 
4 (p 3 -series and interchange the sums to get 

(„-2fe „2fc-|-2.„2v n ( 2m-4^.„2\ (M+4-t.^2\ 

E V9 ;9 i9 jj 2iY^/_.|yl9 j 9 jn-t 19 ,9 )t 

(n'2 n2 n-2m „-2n. ri"2\ A 2-J^ 


(B.IO) ^ (9^ 9"^™, 9"^’^; q^)3 




l^q2m+4. q2-^^_^ (g2;g2)^ 


X (g-2*,g-2—2"+2t-2.g2^^.(^ _ ^2m+2„+2-4i)^2(‘)-(4^+4)t^ 

which has a well-poised structure. Relabeling t = j + p gives that the inner sum is equal to 


(9-2'"-2-2), (g-4^-2;g2)^. 


(^q4e-2m-2n+2.^q2'^^ (^q2.q2'^. 


2j . ^ 2 \ f —2m—2n— 


(9-"^;9^),(9 


2+2j;g2)^.(l _ ^2m-|-2n-|-2-44^^(4^-|-6)4 


^ 3 ( —A£ — 2-\-2j „2j—m—n+1 „29 —m—n+1 2m—2n—2+4?. „2'\ 

W_ -iji _ ^ y _5_y_ ;y Jp (4£-2.7+2 )p 

/—9 ,.AP — 9m—9rj4-94-9.9') — m —n — 1 -.9^—m—r> — 1. .9\ ^ 

p =0 


^g 2 g4£— 2 m — 2 n+ 2 + 2 j q2j — m—n—l q2j—m—n — l.q2'J 


Multiplying the inner sum with ^^^ 2 j- 2 n ’g 2 j- 2 n.’g 2 ) , we can rewrite the sum as a very-well-poised gi/Js-series 

( „ 47 — 2 m — 2 n —2 „ 27 —m—n +1 „27 —m—n +1 2 —4^+2? ^2j — 2n ^2j — 2m 

Q iQ iQ-’ .2 4^-2i+2 

^ 2 j—m—n —1 ^2j—m—n — l f^4i—2m—2n-\-2j-\-2 ^2j — 2m f-,2j — 2n ’ ^ ^ 

y 5 y 5 y ? y ? y 

This very-well-poised series can be evaluated by using [TUI as 

/„4j — 2 m— 2 n „4^— 2 n-l- 2 . „ 2 '\ 

\H 5 y 7 y M—i 

(^4£-2m-2n-h2i-h2^ ^2i-2n. ' 
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Straightforward calculations show that (IB.101) is given by 


(9 


2m-U. „2 




('^2m+4. q,2^ 


(1-9 


2 m+ 2 n +2 


{q-^ra-2n. ^ 2 )^ (^ 4 ^+ 2 - 2 „. ^ 2 )^ ^ 


-2fc ^2fe+2 g-2^-2.^2^ 

(q,4<!-2m-2n+2.q2)^ (q-2n.g2)^ ^ (^2^ ^2 ^ ^,-2^. ^2 




4=0 


The inner sum can be rewritten as a balanced 3 (/? 2 - Using the g-Saalschiitz transformation [101 (11.12)], 
we find that (IB.101) reduces to 


(B.ll) 


(^2m-4^;g2) 

(q2m+4.g2)^ 


■(1-9' 


2 m+ 2 n+ 2 ^ 


/„ — 2 m — 2 n. „ 2 \ /„4^+2— 2 n, „ 2 \ (r^—2k „4^+4. „ 2 \ 

(9 _ i9 jn (.9 ,q )n (q ,q iQ jk 

(^4<'-2m-2n+2.^2)^ 


(^q-2u.q2) 


(92,94^-2'=+2;g2),- 


Finally, simplifying the g-Pochhammer symbols we prove that (IB.Ill) is equal to the right hand side of 
Lemma 16.61 □ 

B.3. Proof of (Tra . Here we verify that (17.31) is valid entry-wise. It follows from Theorem 14. 131 Propo¬ 
sitions [Ol and [OH and the fact that J\fi{z) = that the (m, n)-th entry of the left hand side 

of (17.31) is given by 


(B.12) 


i—m—2n-\-2k-\-3 


/c =0 


(l-9^”)( 


c 


1.-1 71-1 


2 ^—£-\-m-\-n—2k 


-S^+k.e-m-S^+k.m-eJ (1 _ g2^2 


-E' 

fc =0 


m— 2 n+ 2 fe+l 


c 


, 2 / 

• ^ -\-k,i—m — ^-\-k,m—£ 


2 ^^2 ^ 2 n+ 2 ^^—£+m+n— 2 A: 


(l_g2)2(l_^2) 


On the other hand, the (m, n)-th entry of the right hand side of (rr^ is given by 


(B.13) 


E 


^2£—m —2n+4fc+l (2 


(l-g-V) 


— 2 .,^ 2 '^^—£+m+n— 2 A; 


— ^-\-k,£—m — ^-\-k,m—£ 


{l-q^)2{l-z2) 


If we multiply (IB.121) and (IB.131) by (1 — g^)^(l — z^), they become Laurent series in z^. By equating 
the coefficients of z^, the proof of (ESI) boils down to prove the following equality: 


m—2n+2fe+3 


fl _ ^ " 2 ^’^ 'l _ „ 


m —2n+2fc+3 f ^ 22 ’ 


— -\-k-\-l,£—m — ^■-\-k-\-l,7n—i 


m+2A;+3 


_ ^£—m — 2n-\-4k-\-3 


C 


+ 9 

■ ^—m—■j+/c+l,m—^ J 


2 2 ’■*' 

^ — ^-\-k,£—m—^-\-k,'m—£ 


m—2n+4fe+5 


c 


— m — ■^+fc,m—. 




The last equation is proven by using (IA.3I) and performing some simple manipulation of the g-binomial 
coefficients. 


B.4. Proof of (17.5|) . The expressions for /Ci and IC 2 can be obtained using the inverse of L given in 
Theorem 14.151 so that JCi and IC 2 are uniquely determined. So it suffices to check ()7.5I) . 

In order to prove EH) we need to distinguish between the cases i = 1 and f = 2 , since there is no 
symmetry between the two cases. 

The case f = 1 of (17.51) is L*^{z)Mi{z) = /Ci(z)L*(gz). Consider the (m,n)-entry of ICi{z)L*{qz) — 
t{z)Mi{zy, 

(B.14) ^l{z:)7n,mL{qz)n^m 4^1 1 (z). 

The convention is that matrices with negative labels are zero. In case m > n, (IB.14|) is zero, since L is 
lower triangular. In case m = n, we see that ICi{z)m,m = ■Mi{z)m,m implies that (IB.141) is zero as well. 
Note that this also covers the case n = 0. Assume 0 < m < n, and multiplying (IB.14[) by (1 —g^)^(l —z^), 
and using the explicit expressions from Theorems 14.151 and 17.11 and canceling common factors, we need 
to show that the following expression equals zero 

q-m^l _ c„_™(gz; g“|g 2 ) 

+ q-^z (1 - g2-+2™+2) C„_,„_i(z; g“|g2) - g-’^(l - g^+^-z^) C„_™(z; g“|g2). 
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Now we can use the Laurent expansion of (14.101) to rewrite the expression as 


\l-q 


2 + 2 m^ 2 - 


(„2m+2. „2\ („2m+2.„2\ 

E [Q tQ )k\Q ,q )n-m-k f _ s^n-m- 2 k 

k =0 


+ g“” (1 - q 


{q^\q^)k{q^]q^)n-m-k 

n-m-l , 2m+2. „2\ /„2m+2.„2'i 
2ra+2m+2', W_it/ )k\q ,q )n-m-l-k _n-m-2k 


k =0 


{q^] q^)k{q^] q^)n-m-i-k 


l'„2m+2. „2'i |'„2m+2.„2'i 

„2+2n.^2\ yq ly )k\q jy Jn-m-k n-m-2k 

It is a straightforward calculation to show that the coefficient of g |_x^ q, • • • , n — m} equals 

zero. This proves the case i = 1 of (EH). 

To prove the case z = 2 of EH), we consider the (m,n)-th entry of IC 2 {z)L*{qz) — U{z)M. 2 {z), which 
is slightly more complicated than the corresponding case for z = 1 

/C2( —1 ~l~ ^2 ( ~\~ ^2 ('^)7Tl,m+171 ,7?1+1 

(B.15) 

L(z)7i-|-i^j7ijVl2(2^)n+l,n L{^Z^ri^m-^2{,z')n^n- 


If TO > n +1 all terms vanish in (IB.151) . because of the lower triangularity of L. In case to = rz +1, (IB.151) 
reduces to K 2 {z)n+i,n — ■M. 2 iz)n+i,n, which is indeed zero. Suppose m < n, we expand (IB.151) under 
the convention that continuous g-ultraspherical polynomials with negative degree are zero. Expand the 
continuous g-ultraspherical polynomials in z in (IB.151) and take out the term 

( 1 _^ 2 ) (l_^ 2)2 

By a long but straightforward calculation we can show that the coefficient of ^ g {—1,0,..., n— 

TO + 1} equals zero. This proves the case z = 2. 
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